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ABSTRACT 


Highly-uniform, hollow glass spheres (shells), which are used for 
inertial confinement fusion targets, are formed from metal -organic gel 
powder feedstock in a vertical furnace. The modeling of this gel-to-shell 
transformation has consisted of three phases: a study of the effect of 

pyrolysis on gel morphology and thermochemistry; development and testing 
of a furnace-to-gel heat transfer model; and development and preliminary 
testing of a model describing the gravity-driven degradation of shell 
concentricity as a function of shell characteristics and time. 

As a result of the rapid pyrolysis caused by the furnace, the gel is 
transformed to a shell in five distinct stages: a) surface closure of the 
porous gel, b) generation of a closed-cell foam structure in the gel, c) 
spheridization of the gel and further expansion of the foam, d) coalescence 
of the closed-cell foam to a single-void shell, and e) fining of the glass 
shell. Initial foaming is driven by the pyrolysis of the residual organics, 
such as the metal alkoxides and alcohol, and by the decomposition of 
silanols and the alkali bicarbonates in the gel. Further foam expansion 
IS driven by the oxidation of the elemental carbon produced in the initial 
pyrolysis and by decomposition of the alkali carbonates in the gel. The 
sole shell-forming gases were identified from residual gas analyses as 
CO2, H2O, O2 and N2. The O2 and N2 from the furnace ambience permeate 
into the shell during its formation. 

The heat transfer from the furnace to the falling gel particle was 
modeled to determine the effective heating rate of the gel. The model 
predicts the temperature history for a particle as a function of mass, 
dimensions, specific heat, and absorptance as well as furnace temperature 
profile and thermal conductivity of the furnace gas. This model has been 
experimentally verified. 

In the third phase, a model was developed that predicts the gravity- 
induced degradation of shell concentricity in falling molten shells as 
a function of shell characteristics and time. The model predicts 
negligible rates of decentering for glass shells smaller than Z mn 
diameter. The validity of the model and initial efforts for testing 
the model are discussed. 



1 . INTRODUCTION 


Targets currently used in inertial confinement fusion (ICF) include 
spherical glass shells filled with deuterium-tritium fuel. Currently these 
shells have diameters from 100 to 500 ym, and wall thickness from 0.5 to 
10 ym. To meet the experimental requirements the non-concentricityf 

inner and outer surfaces of the shells must be less than 5%, asphericity less 
than 1% and surface irregularities less than 0.5 ym. These specifications 
are becoming more stringent as the available laser power increases and laser 
fusion experiments become more sopnisticated. In addition, the shells are 
required to have high strength, good chemical durability and be permeable 
to various diagnostic and fuel gases. Meeting these requirements becomes 
particularly difficult for the next generation of targets whose diameters 
are greater than 1mm, the effective size limit for current manufacturing 
methods. For such large glass shells, the force of gravity and the aero- 
dynamic forces, acting upon falling molten shells, will increasingly manifest 
themselves, degrading the concentricity and sphericity of the shells, and by 
virtue of the relatively large terminal velocity of the gel, limit the glass 
fining time interval during wnich a falling particle is exposed to the heat 
of the furnace. 

These shell? are currently manufactured by several methods [1,2,3]. One 
of the more promising and versatile methods involves the use of metal-organic 
powder which is fed into an electrically-heated vertical tube furnace whereby 
glass shells are formed. Metal -organic gels as glass precursors permit 
the formation of hollow glass spheres with a variety of glass compositions 
[4,5.6]. 

In view of the limitations on shell size and quality that are imposed 
by the gravitation and aerodynamic forces on the current methods, the 
manufacture of shells in the near-weightless environment of space seems to 
be a viable alternative to earth-based manufacture of large-diameter 
highly-uniform shells. A necessary preliminary to any glass manufacturing 
experiments in space is a better understanding of the gel -to-gl ass-sphere 

^Non-concentricity is defined by the ratio of Ad/d where d is the wall 
thickness. Asphericity is defined by the ratio AR/R, that is the 
difference between the maximum and minimum external radii of the oblate, 
divided by the minimum radius. Surface irregularities are surface 
features, debris, or reaction products deposited on the surface. 
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trans-f^ormation, and an assessment of the effects of gravity on these 
processes and upon the final molten glass shell. Such an understanding 
could be attained by a thorough investigation and modelling of the mechanisms 
tnat comprise the gel-to-glass shell transformation that occurs during the 
terrestrial manufacture of shells. 

Our ground-based research has consisted of three phases: 1) determining 
the gel morphology and thermochemistry, 2) modelling furnace-to-gel heat 
transfer and 3) modelling the effects of gravity on the concentricity 
of the final glass shell. The results of these efforts are presented herein. 



2. PYROLYSIS STUDY OF GELS 


Since the transformation of the gel to glass involves both physical 
and chemical processes a study of this transformation involves both a study of 
the gel morphology and gel thermochemistry. Because the transformation is 
rapid and the furnace temperature high, monitoring in real-time is difficult. 
Consequently, we have performed experiments in which we examined stages 
of the process and typically at much lower heating rates than would occur 
in the vertical furnace. It should be noted that at the lower heating rates 
the gel responds more uniformly to temperature changes and internal stresses 
and the resulting volatiles have time to escape from the gel bulk by 
diffusion. In contrast to this, the rapid heating rates achieved in the 
vertical fi'rnace can be expected to trap more volatiles. This difference 
in heating rate undoubtedly introduces other effects. 

2.1 Gel Morphology . The changes that occur in the structure of the 
gel upon heating were studied by hot stage microscopy. Sodium silicate 
gels (e.g. 20 mole % Na 20 , 80 mole % Si 02 ) were examined extensively because 
they were simple in composition an-' were expected to exhioit most of the 
physical and chemical changes of other glasses of interest. In order to 
compare the sheil-blowing behavior of metal-organic gels to gels having little 
or no organic constituents, two other gel samples were examined. The first 
of these was a metal -organic gel, prepared from sodium and silicon ethoxides, 
which after initial drying, crushing and sieving was exposed to moist air 
for extended periods. We refer to the gel as-made as "unhydrolyzed" and 
that exposed to moist air as "hydro yzed". The second sample was prepared 
by the evaporation of water from an aqueous (inorganic) sodium silicate 
solution. 

When the aqueous inorganic sodium silicate gel was heated it formed 
spheroids at 150°, and crude, transparent glass shells at 300°C. The 
hydrolyzed metal-organic gel, which like the inorganic silicate sample 
had a very low surface area (< 0.1 m /gm), also produced a foam spheroid 
at relatively low temperatures (250 to 400°C), but a similar crude hollow 
shell did not form until about 800 to 900°C. The un hydrolyzed gel, which 
contained significant amounts of organic species, was water-poor and had 
a high surface area ('^^ 300 m /gm) compared to the other samples, remained 
rigid until about 500°C, but sporadic movement of the gel pafticles testified 
to some gas evolution. Near 500°C, the gel began to expand and to carbonize. 



forming a closed-cell foam, and at 800 to 900°C finally formed a crude 
foam or multi-void spheroid. These observations are summarized in Table 1. 

Other, more complex metal-organic gels behaved similarly when examined 
by hot stage microscopy, that is, the unhydrolyzed gels remained rigid 
and unchanged until 500 to 600°C, at which time the organic species 
comprising the gel carbonized. These unhydrolyzed gels then further 
experienced a volume increase due to foaming, much like pop-«.^.n, ■' ’lowed 
by glass flow at the surface at 600 to 700°C. Melting of t "glass »nd 
oxidation of carbon occurred at 800 to 1000°C. The oxidation of the 
trapped carbon probably provided a significant quantity of CO 2 fo<' driving 
the shell expansion to completion in this temperature regime. The 
corresponding hydrolyzed gels generally formed a foam structure and expanded 
at 400°C; otherwise they behaved similarly. Shell intermediates that were 
recovered from shell -forming experiments have yielded further insight into the 
gel morphology. These intermediates, formed from 20 to 180 pgm cylinders 
of compacted gel powder (gel pellets)[7] in a vertical furnace, ranged 
from partially-expanded gel pellets to black spheroids having an internal 
closed-cell structure to crude, carbon-laden spheres. Scanninc electron 
microscope views of these intermediates. Figure la-ld, show the sequence 
of changes in the external and internal gel structure as the gel pellets 
transform into hollow glassshells. At relatively low particle temperatures 
500°C), before significant surface flow occurred, the gel grains that 
comprise the pellet were expanded by the gaseous pyrolysates generated in 
the gel; a closed-cell foam structure was formed in the gel. Upon heating 
the gel to approximately 600®C, further pellet expansion, sintering of the 
grains and significant surface flow occurred, cor-olidating the pellet surface. 
At higher temperatures, as the material became fluid, the surface tension 
became dominant, causing the spheridization of the gel particle and the 
subsequent collapse of the internal closed-cell foam structure into a crude 
hollow sphere (Fig. Id). This consolidation of the internal closed-cell 
structure within the spheroid is more readily inferred from the decreasing 
number of voids in the spheres in Figure 2. The product of this process 
IS a black spheroid whose thick walls contain bubbles and carbon inclusions. 

The last phase of the process Is further expansion and fining of the hollow 
glass sphere, driven by the oxidation of these carbon inclusions. 

Experimental work in support of heat-transfer modelling (discussed later) 
indicated that the transformation of the gel to a hollow black spheroid 

occurs rapidly(w1t*^'‘n M sec 9 1000-1 300®C) for gel particles as massive 
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Table 1. Gel Morphology From Hot Stage Microscopy 


Temperature for Structure Changes (°C) 


Sodium 
Silicate Gel 

Rigid 

Gel 

Sintering 

Foaming 

Expansion 

Spheroid 

Crude Shell 
Central Void 

Aqueous-Based 

25° 

- 

- 

150° 

300° 

Metal -Organic 
Hydrolyzed 

25° 

- 

150-200° 

240-400° 

800-900° 

Metal -Organic 

25° 

400° 

501 

800-900° 
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Figure 1 b. Gel Sintering and Expansion oi Cell Structure. 
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Figure 1 c. Gel Spheroidization 
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as 20 pgm. However, subsequent fining of the carbon and gas bubble 
inclusions in the glass wall was imcomplete even after an additional five 
seconds residence time in the furnace. 

This morphology was deduced from gel particles that were partially 
transformed by heating in a vertical furnace, as well as by supplementary 
experiments in a box furnace. A similar morphology was evident in metal - 
organic gel that was heated to melting while levitated in an air jet. 

This work was done by Dr. Edwin Ethridge of MSEC using gel samples we 
provided [8]. 

On the basis of the above observations from shell -forming experiments 
in our vertical furnace, hot stage microscopy, and a limited number of 
levitation experiments, a general r.iodel for a transformation of tne gel 
morphology can be visualized; a) closur'e of surface of the porous gel 
structure, such as by an initial sintering, b) low-temperature foam 
generation, c) particle spheridization and foam expansion, d) void 
consolidation (foam coalescence) and further expansion to a single-void 
spheroid, an l e) fir.al fining of glass wall of the shell to remove elemental 
carbon inclusions and gas bubbles. The model is a general one, and the 
characteristic temperatures of each state may vary substantially, depending 
on the composition and heating of the gel. 

2.2 Mass/Dimensional Changes . The change in morphology, that is in 
gel 'vass and dimensions, resulting from pyrolysis of the gel was determined 
quantitzbi vely as a function of temperature in the range of 20 to 600°C. 

Since irregularly-shaped gel particles were not suitable for this purpose, 
the unhydrolyzed gel was crushed into a fine powder and compacted into right 
circular cylinders (gel pellets), for which the diameter, length, and mass 
was acc:rately determined [7]. Thus a batch of sodium silicate gel pellets, 
250 n in diameter and 250 pm long, were fully characterized for dimensions, 
•■'ss, composition of pyrolyzable volatiles and specific heat as a function 
of temperature. Platinum crucibles containing 50 mgm samples of these 
pellets were heated in air at 200°, 400°, 500°, 600°, and 800°C for one 
hour; these pellets were subsequently completely recharacterized to assess 
the charges in the gel. The results are summarized in Table 2. This 
dimensional analysis, and subsequent SEM inspection of the pyrolyzed samples, 
shown in Figures 3a and 3b, corroborated the observed initial expansion 
of gel pellets that were pyrolyzed at 400 to 600°C. Pellets that were 
pyrolyzed at 600° had expanded an average Z0% in diameter and 20 to 50% in 
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Table 2. Effect of Heat Treatment on Gel Pellets 


Diameter 


Sample 

Treatment 

(urn) 

1 

untreated 

261 + 2 

2 

200°C 

246 + 2 

3 

400°C 

254 + 5 

4 

500°C 

267 + 2 

5 

600°C 

301 + 6 

6 

800 °C 

1 


Length 

Mass 

Am,%' 

(pm) 

(pgm) 

255 ^ 5 

19.9 + .9 

N/A 

242 + 2 

17.4 + .4 

12.8 

279 + 5 

15.7 + .6 

21.1 

275 + 4 

15.8 + .4 

20.8 

310 + 12 

14.6 + .4 

26.6 


b) b) 


a) Mass loss. Am. is based on the mass (iri) of the untreated pellets 

b) Gel pellets were completely melted and not recoverable as pellets 
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length depending on the extent of pellet densification during compaction, 
and had experienced a corresponding mass loss of 27%. SEM inspection showed 
the pellet surface had formed a continuous outer surface of rough texture 
due to underlying gas bubbles. Enough surface flow had occurred to 
firmly fuse contiguous pellets. 

The 800®C sample had melted completely to puddles of black, foamed 
"glass" that firmly adhered to the platinum boat and were not recovered 
for a determination of the mass loss. Therefore, the mass loss was determined 
for the whole pel let-to-shell transformation, rather than for just the 
molten "glass" at 800°C, from shells formed at IBOO'C from pellets. These 
data are shown in Table 3. These mass correlations were obtained initially 
from discrete pellet-shell pairs. This entailed a) pellet mass and size 
determination, b) formation of the shell from pellet, one at a time, in 
a 1500° furnace, and c) shell mass and size determination. (Once it was 
determined that the gel pellets did not fragment during heating, these 
data were obtained by statistical batch analyses, such as those in Table 4.) 
The average mass loss of 26.2% observed for these shells however, is similar 
to that found for pellets heated to 600° in a box furnace. This indicates 
that most of the mass lost in the gel-to-shell transformation is lost before 
the gel reaches 600°. This conclusion is corroborated by thermogravimetric 
analyses (TGA) of gel powder and gel pellets, wherein most of the evolved 
gases comprising the lost mass are lost prior to 500-600°C, and an 
essentially stable mass results beyond that, up to 1000°. A typical TGA 
thermogram illustrating this for both unhydrolyzed and hydrolyzed gel is 
shown in Figure 4. When dry gel is preheated by heating at 100°C to a 
constant weight to desorb H 2 O, the mass loss observed for the gel by TGA is 
approximately 15% of the dry gel mass. For the untreated dry gel, such as 
the above pellets, typical mass losses observed by TGA are 20-24% (of the 
dry gel mass) somewhat lower than losses resulting from heating pellets 
at 600°C or from the gel-to-shell transformation. Under certain conditions, 
e.g., low gel particle mass, high furnace temperatures and high shell 
surface area to glass mass, it is possible to loose appreciable amounts 
of the alkali and boron present in the glass shell formed from the pellet. 

2.3 Shells From Pyrolyzed Pellets . To determine whether the heat- 
treated gel pellets still contained sufficient blowing agent to form 
shells, the pellets which were pyrolyzed as described above were dropped 
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Table 3. 

Pellet-to-Shell 

Mass Correlations 


Pellet 

Mass (uqm)®^ 

Shell 

Mass (ygm)^^ 

% Mass Loss 

1 

13.5 

1 

9.2 

31 .8 

2 

12.7 

2 

9.9 

22.0 

3 

12.7 

3 

8.8 

30.7 

4 

13.4 

4 

9.6 

28.4 

5 

13.0 

5 

9.9 

23.8 

6 

14.5 

6 

11.3 

22.1 

7 

11.1 

7 

8.4 

24.3 

8 

12.1 

8 

9.0 

25.6 

9 

9.6 

9 

7.0 

27.1 


a) Mass measurement uncertainty + 5%. 
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TGA THERMOGRAMS OF METAL ORGANIC GEL 
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TEMPERATURE °C 



through the furnace at 1300®C (dry air environment). The shell produced 
from the unheated and the 200®-, 400°-, and 500°-pyrolyzed pellets were 
virtually identical in yield, quality, size, and mass, as shown in Table 4 
and Figure 5. The average diameters of the shells from the 500°- and 600°- 
pyrolyzed gel pellets are larger than those for the other pellet samples, 
and also have a correspondingly wider standard deviation for the average 
diameter of the samples. While these data imply that the shell diameters 
increase with increasing pyrolysis temperature, the data actually reflect 
an increasing frequency of shell multipiets ('\.10%) in the shell sample. 

These multipiets had twice the average shell mass; they were dropped 
from the mass analysis. Ho’'»ever their diameters were only slightly larger 

than the normal population and could not be positively identified, nor 
consequently eliminated from the size analysis on the basis of diameter.. 

They therefore skewed the analysis to a "larger ave»^age. It is presumed 
that these multipiets are the products of collision between gel pellets 
or shells in the furnace. Evidently the gases evolved by the pyrolyzing 
gels effectively disperse a gel charge in the furnace, reducing the chance 
of collision. A significant reduction of these evolved gases by prior 
pyrolysis of the gel results in ineffective particle dispersion and an 
increased collision frequency. The 600°C-pyrolyzed sample sintered and fused 
at contact points. These agglomerates, which were used without separating 
them, expanded at 1300°C, but due to their greater mass, formed only crude 
spheroids. The 800°C sample had formed "glass" puddles and was not recoverabl 
Thus the heat treatment, which stripped the gel of considerable amounts 
of volatiles, did not alter the pellets in any way that affected the shell- 
forming process. Evidently the residual elemental carbon and/or the Na2C02 
in the gel (see below) are the sources of CO 2 sufficient to form the shell. 

2.4 Specific Heat . The specific heat of the pellets, determined 
by differential scanning calorimetry, changed markedly as the organic 
and inorganic species in the gel were removed by pyrolysis. Figure 6 
shows the specific heat (Cp) as a function of temperature for two gel 
samples pyrolyzed at 400°C and 600°C, and for the unpyrolyzed gel. The 
figure also includes the Cp trace for the sapphire (AI 2 O 2 ) calibration 
standard. The unpyrolyzed gel ’n general exhibits a higher Cp due to the 
latent energy due to the silanol condensation, organic pyrolysis and NaHCO^ 
decomposition reactions. But superimposed on this is a strong endotherm 
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Figure 5. Shells From Pyrolyzed Gels. 
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EFFECT OF PYROLYSIS ON THE 
SPECIFIC HEAT OF GEL 



TEMPERATURE 








at approximately 100°C, corresponding to the evaporation of adsorbed H 2 O 
and free ethanol, and exotherms at approximately 400°C and 600°C, 
corresponding to oxidation of evolved orrganic pyrolysates and gel sintering 
reactions. These are totally absent in the Cp trace for the gel pyrolyzed 
at 600°C. The trace is similar to that of the sapphire standard (though 
somewhat endothermic) and to Cp values for glass. Gel pyrolyzed 
at 400°C exhibits the expected features, i.e., below 400°C it exhibits 
Cp values similar to sapphire and the 600°C-pyrolyzed gel, but above 400° 
exhibits the oxidation and sintering exotherm, similar to those of the 
unpyrolyzed gel. The pyrolysis of a gel produces a decrease in the specific 
heat, as measured at a sample temperature of 100°C, from 5.4 j/gm K to 8.4 
j/gm K, due to the elimination of the endothennic pyrolytic reactions. It 
produces an increase in the specific heat at 600° from 0.85 j'/gm K to 1.27 
g/gm K due to the elimination of the exothermic reactions. Thus the specific 
heat of gel samples pyrolyzed at 600° and stripped of all endothermic and 
exothermic reactions ranges from 0.84 (at 100°C) to 1.27 j/gm K (at 600°C). 
This is comparable to that of glass, which ranges from 0.92 (at 100°C) to 
1.30 j/gm K (at 600°C). 

2.5 Gases Evolved by Gel Pyrolysis . Gases are generated in and evolved 
by the gel during pyrolysis, and are the agents that expand the closed-cell 
foam structure that eventually results in a hollow glass shell. As is evident 
from the above mass analyses and shell -forming experiments, most of these 
evolved gases escape from the gel, especially during the early stages of 
the pyrolysis prior to sintering of the porous gel. They presur.iauly comprise 
most of the observed mass loss from the gel. These gases were identified 
by pyrolysis-gas chromatography (PGC), pyrolysis-mass spectrometry (PMS) 
and, indirectly, by differential thermal analysis (DTA). The identities 
of evolved gases obtained by PGC analyses, and the typical relative amounts 
of the evolved gases, are given in Table 5. Typical chromatograms of 
the evolved gases are shown in Figure 7 for sodium silicate gels made from 
aqueous silicate solution and from metal-organic solution. The traces also 
show graphically the difference in the evolved species from metal - 
organic gel that has been maintained dry (unhydrolyzed) and that exposed 
to ambient humidity (hydrolyzed). Thus, the evolved gases from all these 
gels consist typically of CO 2 and H 2 O; those gases evolved from unhydrolyzed 
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b) Gels had been dried to constant weight at 100“ prior to analysis. 

c) Samples 1 and 2 are two separate batches of gel synthesized and stored under dry air. Sample 3 was 
synthesized and stored under dry argon to exclude the C0„ in air. 







metal-orgamc gel also include ethanol and its pyrolysates. 

The source of the watef consists of chemically bound species as well 
as molecular water adsorbed to the highly-active gel surface. Chemically 
bound water is present in the form of metal hydrates, silanol groups and 
metal -bi carbonates (see below). The organic gases and additio’’.! ^^tr.r 
arise from the pyrolysis of residual organic groups, such as silicon 
ethoxide in the gel network: 

0 

^ I , 

-tSi- Of + 

0 

+ CH3CH.OH + CH^ = CH^ + CH3CH2 + CH^ 

These gases serve as initial gel ■''oaming (or expansion) agents even before 
oxidation occurs. Complete oxidation of the ethoxides (and the hydrocarbon 
pyrolysates) would further proauce 2 moles of CO2 and 3 moles of water from 
every equivalent ethoxide. The oxidation of the organic pyrolysates (as 
they are evolved by the gel, not within the gel) is suggested by DTA aiialyses 
such as shown in Figure 8. Analyses using air as the ambient gas yielded 
exotherms at ^450®C and at '>^800°C, whereas with nitrogen endotherms appear 
at these temperatures. These results attest to the evolution of oxidizable 
gases; however since elemental carbon is formed within the gel under these 
conditions, oxygen permeation of the gel must be minimal and any oxidation 
must be occurring external to the gel. 

A correlation of TGA, DTA and PMS analyses of such an unhydrolyzed 
metal-organic gel, as in Figure 9, clearly shows the significant mass loss 
due to an initial loss of H^O, ethanol, and CO^ and a rather sharp mass loss 
at approximately 450° due to the loss of the organic pyrolysates. This 
mass loss at 450° corresponds to the rather rapid formation of a closed-cell 
foam structure in the gel at 400-500°. 

The pristine metal-organic gel, which is produced by a controlled partial 
hydrolysis of the precursor solution (in ethanol) of the appropriate metal - 
alkoxides, contains a significant number of ethoxide moieties in the dried 
gel. The moieties are expected to be water sensitive. While the dried gel 
has been shown to be stable when stored under anhydrous conditions, the gel 
does age when exposed to ambient humidity, with a consequent reduction in the 


0CH,CH, 

I ^ ^ 

0-(Si-0f 
I n 
OH 
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TEMPERATURE 




JURE 






ethoxide content. In fact, removal of the bulk of the organic residues 
can be achieved by controlled hydrolysis of the gel through exposure of 
the gel to a constant high humidity environment. The extent of hydrolysis 
and tne reduction of the organic constituents in the gel are reflected in 
the data of Table 5. While they show a decrease in organic pyrolysates, 
hydrolyzed gels also show an increase in evolved CO^ in pGC analysis. 

Since the gel is pyrolyzed in this analysis in an inert helium atmosphere, 
this suggests that the evolved CO2 is not organic in nature. This 
conclusion is substantiated by the considerable CO2 levels in PGC analyses 
of sodium silicate made from an inorganic, aqueous solution. A source 
of the evolved CO2 may be alkali carbonates or bicarbonates in the gel. 
Adsorbed or trapped air in the porous gel will not supply enough CO2 to 
blow shells since the natural level of CO2 in air is only 0.033 percent. 

A CO2 concentration mechanism is required, and can exist in the form of 
the reactive alkali (sodium) in the gel. Sodium bicarbonate, NaHCO^, and 
sodium carbonate, Na2C02, can form upon exposure of the moisture-sensitive 
metal-organic gel to the humidity and CO2 in the air. These reactions are 
thermally reversible, and should yielded CO2 in two distinct temperature 
regimes: a) below 300°C, by decomposition of NaHCO^ to Na2C0^ and b) 
above QOO'^C, by decomposition of Na2C02 to Na20. When the gel was synthesized 
and handled under argon to prevent the formation of the carbonates due to 
exposure to ambient air, the PGC analysis showed a significant reduction 
of the CO2 levels in the gel pyrolysates (Table 5). 


29 


Iki 



r 


3. NATURE & SOURCE OF SHELL FORMING GASES 


3.1 Residual Gases In Shells . The intermediate stage in the gel- 
i to-glftss shell transformation, a '•losed-cell foam structure, is formed by 

the gases generated and trapped in the pyrolyzed gel. These initial blowing 
gases are H^O, CO 2 and the hydrocarbons methane, ethylene and ethane. The 
final stage, the fined hollow glass shell, is formed by the blowing gases 
^ that remain in the shell as residual gases. These gases were presumed 

to consist primarily of the descendents from the above pyrolysates and 
to originate froin the organic species in the gel. The composition of these 
residual gases in the glass shells was analyzed by gas chromatography; the 
identification of the components was confirmed by mass spectrometry. 

Carbon dioxide is the major gas constituent in the shells prepared 
from metal -organic gels. Oxygen and nitrogen are secondary, but together 
comprise 10 to 40 percent of the residuf. '»s. Other possible permanent 
gases, such as carbon monoxide or riitrc -'S, were not detected by 

either gas chromatography or mass spectr . The relative concentrations 
of these gases in four different shell samples, as determined by gas 
chromatography, are presented in Table 6. These data suggest a correlation 
between residual gas composition and glass composition but the information 
presently available is insufficient to establish this relationship. The 
pressure of the residual gases at room temperature is 0.17 to 0.20 atm 
for three of the samples; this is the expected range for shells initially 
at pressure equilibrium at the manufacturing temperature of 1200-1 500°C 
and subsequently cooled to room temperature. The residual gas pressure in 
the calcium containing glass, however, is higher than would be expected 
from similar calculations based on the ""urnace temperature and glass 
surface tension considerations. This high residual pressure is not an 
artifact, having been observed in several calcium-containing glass 
samples. 

Water was analyzed but does not appear in the gas analysis results 
presented in Table 6. Nevertheless, it is believed to be an important 
species in shell blowing, and, depending on its concentration in the gel 
and the furnace ambience, can influence several glass properties and shell 
characteristics. Indirect evidence for the existence of water vapor in the 
shell during the blowing process is the presence of "weathering" products 
detected on the interior surfaces of most shells by scanning electron 
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Table 6. Residual Gases in Shells of Several Compositions 



Relative 

Concentration 


Shell Glass 

Volume 

%, 

By GC 

Pressure 

Sy sterna) 

CO 2 

“2 

"2 

(atm) 

Na.Si 

(Unhydrolyzed) 

71 

19 

10 

0.19 

Na, Si 

80 

13 

7 

0.175 

(Hydrolyzed) 
Na, K, B, Ca, Si 

63 

16 

21 

0.23 

(Unhydrolyzed) 
Na, K, B, Si 
(Unhydrolyzed) 

88 

8 

4 

0.20 


a) The metal oxide glasses were prepared from dried 
metal -organic gels which were either subsequently 
kept dry (unhydrolyzed) or exposed to 80 % RH for 
about 5 days (hydrolyzed). Shells were blown in 
a dry air environment. 
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microscopy. These features are typically rich in alkali (such as sodium) 
and nearly identical in appearance to those found on the exterior of 
shells which are commonly associated with any alkali-rich glass surface 
exposed to water vapor. Water vapor, as high as 30% of residual gases, has 
indeed been detected as one of the residual gases when the shells were 
prepared in an environment enriched in water vapor (See Table 7) [9]. 

When glass shells are formed in furnace environments having low water- 
vapor content, water initially in the gel may participate in the blowing 
process, but subsequently it either diffuses through the shell walls and/or 
is totally consumed by the glass surface and it thus not found in residual 
gas analysis. This can also explain the low incidence of surface features 
on shells blown in dry air. 

3.2 Source of the Gases . The pyrolysis of metal -organic gels generates 
up to 20% by weight (of the dry gel) of volatiles consisting of CO 2 , H^O, 
and ethanol and its pyrolysates as noted above. All of these are potential 
shell -blowing agents. 

An ample source of CO 2 (and H^O) is the oxidation of the ethoxides, 
trapped ethanol and the pyrolysates of ethanol trapped in the gel. If the 
oxidation of these organic residues in the gel is the prime source of the 
COn blowing agent, then a reduction of these residues should affect shell 
formation. A reduction of organics was achieved by controlled hydrolysis 
of dried gel, which cleaved residual ethoxide groups from silicon and 
liberated ethanol. In addition, for comparison purposes, a gel sample 
was prepared from an aqueous sodium silicate; there were no organic 
materials used in the synthesis of the gel. 

Subsequent shell-forming experiments using unhydrolyzed metal -organic 
gel, hydrolyzed metal -organic gel and inorganic gel produced shells of 
supri singly good quality and yield for all three gels although the shells 
from the unhydrolyzed gel had higher aspect ratios (defined as the ratio 
of the shell diameter to wall thickness, OD/W). Likewise, gel pellets that 
were pyrolyzed at temperatures up to 500°C so that only low levels of 
elemental carbon and sodium carbonate remained, resulted in shell yields 
and quality similar to those produced by unpyrolyzed gel, as discussed 
above. These results are contrary to those expected from a comparison 
simply of the potential organic volatiles of these gels, as summarized 
in Table 5, and suggest that oxidation of organic species may not be the 
dominant source of CO 2 for forming shells. 
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Table 7. 


Residual Gases in Shells Exposed to Water Vapor 


During Shell Blowing 


' Relative Concentration, 

Volume %, by GC 

Shell Glass 


Sample 

System 

Furnace Gas 

CO2 

“2 

"2 

H^O 

1 

Na, K, B, Ca, Si 

Dry Air 

57.9 

18.6 

23.5 

0 

2 

Na, K, B, Ca, Si 

Water Vapor, 
Air 

21.3 

14.7 

30.9 

33.1 

3 

Na, K, B, Si 

Dry Air 

88.2 

7.5 

4.3 

0 

4 

Repass of 3 through 
furnace 

Water Vapor, 
Air 

44.3 

17.7 

7.2 

30.7 

5 

Na, K, B, Si 

Water Vapor, 
Air 

35.8 

17.0 

15.7 

31.4 

6 

Repass of 5 through 

Dry Air 

54.8 

33.7 

11.5 

0 


furnace 
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In addition, the PGC analyses of gels, in which the gel is pyrolyzed 
in a helium flow, show substantial levels of CO2 in the evolved gases and 
suggest a CO2 source in gels that is inorganic, e.g., alkali carbonates. 
Minimizing the carbonate formation by synthesizing, handling and hydrolyzing 
the gel under argon, free of CO2, did decrease the available CO2, but 
did not eliminate it. Such gel, i.e., with minimized carbonates and 
residual organic content, produced shells in >ield, aspect ratio, and 
quality cc .oarable to those shells from the three gels above. The results 
show essentially that minimizing the carbonates and residual organics 
in the gel has no significant effect on CO2 levels in the shells. 

Other experiments performed ii, our laboratory have further substantiated 
our conclusion that residual levels of elemental carbon, formed by 
pyrolysis of the residual organics, and alkali carbonates generate 
sufficient CO2 to blow good shells [9]. Only gel pyrolyzed at 900'^C, in air, 
suffers enough reduction in the elemental carbon and the alkali carbonate 
content to result in a significant decrease in the aspect ratio of the 
shell s. 

The source of the O2 and N2 in the residual blowing gases of the shells 
is probably the ambient furnace air, which permeates through the molten 
shell wall, driven by the partial pressure differential between the furnace 
ambience and the shell interior. Significantly, the oxygen is found in 
anamolously high levels compared to nitrogen; whereas in air O2/N2 = 0.27, 
in shells it is often O2/N2 > 1. This suggests that O2 and N2 enter the 
shells by permeating the molten shell wall, but at unequal rates, and are not 
the result of air encapsulation caused by gel sintering. A further though 
less convincing argument against encapsulation of ambient gases, such as 
air, in the gel during shell blowing is the observation that shells from 
gel made and stored under argon, but blown in an air-filled furnace, contained 
air, but not argon. The higher than expected O2 levels are not the result 
of H20 thermolysis, as an increase of the water vapor pressure in the furnace 
decreases, rather than increases, the O2/N2 ratio. 

A further confirmation that permeation of ambient furnace gases can 
occur during shell formation was shown in a shel 1 -blowing experimerit using 
a partial backfill of argon in the furnace. A small amount of argon ('^5%) 
was found in the shells. Shells formed in air, when repassed through the 
furnace partially filled with argon, had similar argon levels as the above 
shells made in argon, and in this case argon could only be present by 
diffusion through the shell wall. 
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4. FURNACE-TO-GEL HEAT TRANSFER MODEL 


The gel powder, falling at terminal velocity through the hot furnace, 
is transformed into glass shells provided that sufficient heat is transferred 
to the particle during the time it resides in the furnace. If the powder 
particles are too massive and the total residence time of the pov^der is 
consequently too short, a product of diminished quality with respec, to 
sphericity, uniformity, and gas bubble inclusions results. With our 3.9 m 
1500*0 furnace, the particle mass limit is of the order of 50 ugm. The 
process is the»"eby currently limited to producing high quality shells with 
diameters that are less than 1 mm. 

This limitation could be the result of several Tactors: a) very slow heat 

transfer from the furnace wall to the gel (poor heat conduction by the furnace 
gas), b) poor absorptior- of the radiant energy by the gel (or high 
reflectivity) , c) slow rates for the chemical reactions that produce blowing 
gas and form molecular bonds in the glass, and d) slow rates of mass flow 
in the gel/glass (high viscosity and slow response to the stress created 
by the blowing gases). These will result in a very slow expansion of the 
particle, a correspondingly slow increase in particle ouoyancy in the furnace 
gas, and relatively short an^ insufficient particle residence times '.n the 
furnace. 

ihese possible limiting factors for shell product‘s on cannot be readily 
investigated and corroborated enipiriu?lly in t>oi ^ ,e in the furnace, due 
to the hostile environment of the furnace, the small mass and dimensions 
of the sample, and the dynamic nature and short duration of tne whole process. 
One of these factors, the possible limitation on shell size imposed by heat 
transfer to the gel, was approached with the aid of a forced-convection heat 
transfe'^ model. This model was developed to calculate the temperature history 
of a gel particle falling through the furnace as a function of gel particle 
parameters and furnace parameters. Thus, for a given set of gel and furnace 
parameters, the model can predict whether the falling gel particle can reach 
the appropriate temperature for a pre-designated transformation, the rate 
at which that temperature is reached, and the residence time in the furnace. 

It does not address the question of v^hether a particular gel particle will 
form a hollow glass shell. 
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4.1 Mathematical Model . The heat transfer model has been written 
into a computer code, termed GELSHEL. It is a composite of two sub-models: 
a) a particle displacement model that calculates, as a function of time, 
the position of a gel particle in a drop-tower furnace, and b) a forced- 
convection heat transfer model that describes the rate of heat transfer 
from the furnace to the gel particle. The GELSHEL code thereby calculates 
the temperature of a gel particle (and the resulting shell) as a function 
of its drop time and vertical displacement in the furnace. 

The pai^ticle displacement model is based on Newton's Second Law 
and Stokes' Law . The velocity of the particle is determined as function 
of time from Newton's Second Law. Its displacement in the furnace is 
determined by integration of the velocity with respect to time. The 
displacement model that describes the particle position within the furnace 
as a function of the time of fall is summarized by: 


dV 

dt 


g (1 - ^) - 12.5 


^f A 
Dm 


V 


dt 


V 


where 

V = particle velocity 

S = particle displacement (distance of fall) 
t = time 

g = acceleration of gravity 
= density of the furnace gas 
= absolute viscosity of the furnace gas 
p = particle density 
m = particle mass 
A = frontal area of the particle 

D = characteristic dimension in Reynolds' number; for the shell 
D is the shell diameter; for the gel, D is the arithmetic 
average of a gel pellet's diameter and length 

The heat transfer model describes the rate of heating of the particle 
in the furnace as the sum of both convective and conductive heating by 
the hot furnace gas and radiative heating from the furnace tube walls. 
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The model is based on the fundamental heat transfer relationship (see 
Appendix I). 





- T) + 


aEA$ 

mC 




where 

m = particle mass 
C = specific heat of the particle 
T = particle temperature 
= furnace temperature 
t = time 

A^ = particle surface area exposed to the thermal environment 
h = convective heat transfer coefficient 
a = Stefan-Boltzmann constant 
E = absorptance of the particle 

The particle, as it drops through the furnace, changes in mass, size, 
geometry and specific heat. This evolution causes substantial changes in 
the heat transfer rate and particle velocity, as a consequence of the changing 
geometry, surface area and mass, and thus particle buoyancy and drag. To 
model these changes, this complex particle history is conceived as occurring in 
three discrete regimes on the basis of particle geometry, wherein the particle 
geometry is determined by the particle temperature. Particle characteristics 
are held constant within each regime, but can be changed from one regime to 
the next in a step-wise fashion. The temperature bounds (temperature 
switches) of the regimes and the particle characteristics within each 
regime are the variables that must be defined by the user. These regimes 
in their chror.wlogical order are: 

1) cylindrical gel pellet regime ; mass, dimensions, and specific 
heat may be changed at the end of this regime, but not the pellet 
geometry, which is a right circular cylinder. 

2) pellet sintering regime ; at the end of this regime, mass and 
geometry are changed (from a cylinder to a sphere), and specific 
heat and absorptance may be changed within this regime. 

shell regime ; absorptance may be changed within this regime, but 
mass, geometry and specific heat are constant for the glass shell. 
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The modelling approach is based on four simplifying assumptions: 1) gel 

pellet characteristics change ' a step-wise and not a continuous fashion, 

2) the geometry of the gel pan.icle is a right cylinder, 3) the geometry 
change from pellet to shell is instantaneous and 4) the heat transfer 
from the surface to the bulk, of both pellec and shell, is instantaneous, 
so there is no temperature gradient within the particle. 

The development of this model is described in detail in Appendix I, 
and the description of the 6ELSHEL code is given in Appendix II. 

4.2 Model Predictions and Experimental Observations . Those parameters 
that can be controlled or determined empirically, or which are expected 
to have a profound influence on the history of the particle, were selectively 
varied in the code calculations. The particle temperatures, however, which 
served as gel regime switch points were determined from TGA and box-furnace 
test data, and were used as such, since variations on the order of 100°C 
did not have a significant effect on the calculated gel history. 

Calculations using the GELSHEL code indicate the following: 

1) An increase in mass significantly decreases the particle heating 
rate. 

2) An increase in particle length, that is commensurate with an increase 
in particle mass, minimizes the terminal velocity, maximizes the 
heating and results in faster particle heating in comparison to a 
comparable increase in particle diameter. 

3) Heat conduction from the furnace ga^^ is the primary means of heat 
transfer, especially at temperatures lower than 1000°C. Not only 
does the heat conductivity of the gas determine its efficiency as 
a heat exchanger, but the gas viscosity and density affect the 
buoyancy and drag on the falling particle as well. 

4) The initial specific heat of the gel has only a modest effect on 
the pellet temperature history. 

GELSHEL code calculations indicate that an increase in the gei pellet 
mass results in a significant increase in the droptime and distance required 
by the gel geometry changes that accompany the mass increase. If the pellet 
diameter is doubled, commensurate with a four- fold increase in mass, and the 
length is held constant, the drop distance required by the pellet to reach 
1000° is increased six-fold. !f the pellet length is quadrupled to inci'ease 
the mass four-fold, and the diameter is held constant, the drop distance is 
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only doubled. This behavior follows from the following considerations. 
Since the rate of particle heating is given by: 


dt 


Agh (T^ - T) 
itiir 


Ag aE - T^) 

mCT 



+ 7iDL 


m 


= Vp = 


irD^Lp 


then the change in dT/dt for a change in L is given by: 


3 (^\ -2 a 

3L Vdty " 

h (T^ -T) + o E (t/ - 
Similarly for a change in D 


where a = 


pu 


-L lAa_ 

DC Vdt/ ■ p2 

Thus, as expected, any increase in dimensions in the pellet with constant 
density, will result in a decrease in the rate of temperature rise. Further, 
the rate of change in -^will be the same for any change in D or L if the 
ratio of the pellet dimensions ^ = / 2 
That is. 



Since the ratio of dimensions for the pellets that can be manufactured with 
our pellet press lies in the range 
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0.5 < ° <2 

the above relationship and its inverse, i.e., ^ ^ 2 , and 

IL„ 1< ll_ (^\\ 

|3L \dt/ I I9D Vdt/1 , both apply for the range of pellets in use. 
Further, the particle velocity is also affected by the pellet dimensions 
and mass because the buoyancy and drag terms in the particle acceleration 
equations are geometry-dependent. With proper definitions for the buoyancy 
and drag terms, the acceleration for the pellet is give by 


dt 


= g ( 1 - 


D^L 

4m 


12.5 


Pf LV 
m 


Partial differentiation with respect to 0 and to L yields 



Thus, for a given increase in pellet mass with a commensurate increase in 
diameter, the decrease in pellet acceleration is a function of the diameter 
and mass. For a commensurate increase in pel let length, the decrease in pellet 
acceleration is solely a function of the mass. Since the pellet mass 
and dimensions are interrelated variables, the above considerations 
significantly complicate an analysis of the effects of pellet mass on the 
thermal history of the pellet. The parameter of pellet mass in conjunction 
with pellet dimensions is nevertheless a sensitive probe for testing the 
heat transfer model. 

Accordingly, a pellet series was fabricated with a nominal mass of 20, 

3 

40, 80, 160 and 200 pgm and a density of 1.4 gm/cm . Because the pellet 
diameter is fixed by the press die dimension, it could not be continuously 
varied but was limited to two values - 250 ym and 500 ym. The pellet 
length was variable but in a relatively narrow range from 0.5 0 < L<1.5 D. 

Furnace experiments were conducted with several furnace profiles: 
a) 1500°C, all units; b) 1300°C, all units; c) zoned, 1300°C top unit and 
remainder at 800°C; d) zoned, 1200°C top unit and remainder at 800°C; and 
e) zoned, 1000°C top unit and remainder at 800°C. 

The results are best exemplified by a summary of the observed and 


40 



calculated results for experiments using furnace profile b, and are shown 
in Figure 10. The drop distances for the pellets were visually estimated 
at the point that the particle incandescence became indistinguishable from 
the ambient furnace radiation. Because of the speed of particle heating, 
the small size of the particle and the subjective nature of this determination 
of particle temperature equivalence, the particle temperature at this point 
is most likely in the range of 1000°-! 300® C. Therefore, Figure 10 includes 
the calculated times and distances (by the GELSHEL code) that are required 
for the particle to reach the temperatures of 1000° and 1200°C respectively. 

The relatively minor differences between drop times and the corresponding 
particle drop distances for 1000°C and 1200°C illustrate the difficulty 
in estimating the point at which the particle is approximately equal to the 
furnace temperature in the visual tracking technique. Nevertheless, the 
calculated distances required for each pel'et mass to reach 1200°C aro in 
remarkably good agreement with the experimental observations. 

From model calculations, the specific heat of the gel was expected to 
exert some influence, albeit a moderate one compared to gel mass and furnace 
gas, on the thermal history of the gel. The GELSHEL model predicts that 
a decrease in the specific heat sfiould increase the pellet heating rate and 
result in a shorter drop distance (into the furnace) to bring the particle 
to its melting temperature. The thermal histories for gel pellet samples 
were calculated by the GELSHEL code, using the appropriate pellet mass and 
dimensions, for three different furnace temperature profiles: a) all furnace 

units at 800°C, b) top unit at 1000°C, and the remaining 5 at 800°C, c) top 
unit at 1200°C, and the remaining 5 at 800°C. In the profiles b and c, 
the hot zone (1000° and 1200° respectively) at the top of the furnace was 
50 to 60 cm long; at 100 cm depth, the temperature had reached the 800°C 
baseline. 

The effect of the lower specific heat of the heat-treated pellets is a 
faster heating rate. However, this effect is felt for only a short temperature 
range, since the specific heat rapidly decreases to that of glass as the 
gel is heated to 600°C and higher. This effect decreases with increasing 
furnace temperatures. The specific heat affects primarily the early portion 
of the pellet temperature history and to a lesser degree the maximum particle 
temperature, the total droptime, and distance required to attain it. 

These data are summarized in Table 8. These calculations are in good 
agreement with the drop distance of 30-50 cm observed for these in a 
furnace with temperature profiles b and c. 



PELLET DROP DISTANCE AS A FUNCTION 
OF PELLET MASS TO REACH 
FURNACE TEMPERATURE 
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Table 8. Effect of Specific Heat on Pellet Heating 


Gel 

Treatment 

Specific 


Peak 

Peak 

Particle Temp. 

(3) 

Heat^is 

Furnace/^v 

Profile'^'' 

Fu rnace 
Temp. °C 

Temp. °C 

Time, s 

Distance, cm 

Untreated 

5.44-0.85 

a 

1045 

991 

0.52 

43.8 

400“C/1 hr. 

1.31-0.70 

a 

1045 

1024 

0.46 

33.1 

500°C/1 hr. 

1.19-0.25 

a 

1045 

1025 

0.45 

32.5 

Untreated 

5.44-0.85 

b 

1202 

1161 

0.43 

27.7 

400°C/1 hr. 

1.31-0.70 

b 

1202 

1179 

0.42 

27.7 

500°C/1 hr. 

1.19-0.25 

b 

1202 

1180 

0.42 

28.2 


(1) Range of values, measured between 100°C-600°C 

(2) a) top unit set at 1000°C, lower 5 at 800°C 
b) top unit set at 1200°C, lower 5 at 800°C 

(3) Drop time (secs) and drop distance (cm) required by 
particle to reach the stated temperature based on 
calculations using GELSHEL code. 
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The third probe for testing the thermal model is the thermal conductivity 
of the furnace gas. The thermal conductivity, gas density and gas viscosity 
are expected to determine the rate at which thermal energy is transferred from 
I the furnace wall to the pellet, as well as the velocity of the pellet 

through the heated zone. The 6ELSHEL code predicted that only helium was 
significantly more effective for heat transfer, and water vapor only 
moderately better than dry air. Because developing and maintaining a 
helium atmosphere in the open process tube of the furnace would require a 
fast gas flow, thus perturbing the particle trajectory, water vapor was chosen 
for the test. The experiments were performed with the heavier pellets to 
minimize any perturbation of the pellet trajectory by the updraft of water 
vapor and thus to generate a more valid comparison with the dry-air control 
experiment. These results are presente'^ in Table 9. The effect of the 
higher heat conductivity of water vapor is to shorten the drop distance 
required to heat the gel, and is seen in both the GELSHEL calculations and 
the experimental observations. Further, the dry air and the water vapor 
experiment observations and corresponding calculations show a remarkable 
agreement between the model and the real particle behavior. In these 
experiments, any great disparities between the calculated and observed 
results probably reflect the inaccuracy inherent in the visual estimates 
of distance and particle temperature, especially in those cases wherein the 
particles plummet deep into the furnace before reaching the furnace 
temperature. 

In all furnace experiments, to date, the pellets reached incandescence 
very rapidly corroborating the GELSHEL calculations, yet the final product 
from the more massive pellets did not have the surface characteristics 
expected for such a hot thermal history. The physical appearance of some 
of the larger (180 ygm) pellets, dropped through a zoned furnace, was 
equivalent to that of reference samples equilibrated at 500-80o’C in a box 
furnace. Further experiments (in a box furnace and vertical furnace) wherein 
pellets were heated to 1000°C and 1200°C for short periods (< 5 secs), 
indicated that while the pellets reach incandescence rapidly, the mass flow 
is negligible. There are two explanations for these observations: 

1) While the thermal response of the pellet to the furnace 
environment is rapid, the viscous response of the gel to 
its internal stresses is much slower in comparison. 
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2) The pellet, because of limited heat conductivity, may develop 
a significant temperature gradient resulting in an incandescent 
surface and a substantially cooler bulk. 

This latter point questions the validity of a basic assumption in the 

GELSrIEL model. Namely, because cf the small pellet mass, the heat transfer 

in the pellet bulk is assumed to be instantaneous; no significant temperature 

gradient develops in the pellet. The radial temperature gradient was 

ther:fore calculated for a typical pellet with a 20 ygm mass and a 300 ym 

diameter, having the thermal properties of glass (specific heat of O.S-1.3 

_2 

j/cm K; heat conductivity of 10 w/cmK) as discussed in Appendix III. 

The heat flux to the pellet, due to convection and radiation, was calculated 
from the heat transfer equation. For a pellet injected into a furnace having 
ideally a constant temperature profile of 1500°C, the maximum radial 
temperature gradient of 75°C developed within 4 msec droptime. This 
gradient rapidly decreased beyond that time, to 50°C after 20 msec, 
corresponding to a pellet displacement of only 2 mm into the 1500°C furnace. 
Since in reality the mouth of a 1500°C furnace is only 750°C, a pellet that 
has descended 2 mm into the furnace sees a 750°C environment. A recalculation 
of the radial temperature gradient for such conditions resulted in a 22® 
gradient that decreased to 12 after 20 msec (see Appendix III). Similarly, 
a 180 ugm pellet develops a maximum radial temperature gradient of 105°C 
within 20 msec after injection into a 1500° environment, and 25°C within 
20 msec after injection into a 750° environment. These radial temperature 
gradients are indeed not significant in view of the fact that they 
maximize within a very short time after injection, and rapidly decrease 
to lesser values. 

To differentiate empirically between a temperature gradient or viscous 
inertia in the pellet requires a model material with a specific heat similar 
to glass, but with a high surface tension, a sharp melting pcint, and very 
low viscosity (for the molten mass) to drive the rapid spheridization of the 
molten mass. Thus, such a material would heat like glass but spheridize 
rapidly at its melting point. Except for their high heat of fusion, inorganic 
salts, with a melt viscosity near 1 centipoise, seemed suitable. Accordingly, 
NaCl and Na2S0^ crystals of 212-250 urn sieve range and 20-25 ugm mass, were 
injected in the furnace in the usual manner. With the top unit of the furnace 
at 1000°C and the remaining units at 750°C, NaCl crystals (m.p. 801°C, 
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AH^ = 6.73 + 0.04 kcdl/mole) became incandescent within 30 cm drop; the 
product consisted largely of solid spheres, with some unaffected crystals. 

With the top unit at 900°C, NaCl did not spheridize even though incandescence 
was reached and maintained for 1 m drop. This failure to spheridize, even 
with exposure to at least 65 cm of travel through a furnace zone above 
800°C, reflects the high heat of fusion of NaCi (6.73 kcal/mole) and the 
short residence time, rather than a significantly large AT in the particle. 

Had the product experienced a large temperature gradient within the particle, 
a substantial number cf partially-melted crystals should have resulted, 
whereas the products were sharply divided into spheres and unaltered crystals. 
Similar results were obtained for Na2S0^ (m.p. 884°C) fired at 1200° and 
1000°C. 

The qualitative agreement of the experimental results with the 
calculations from the model verifies that the model can satisfactorily 
describe the temperature history oi the gel in the furnace. Thus, a gel 
particle with a mass of 20 ygm and a diameter and length of 250 ym is 
heated rapidly to the temperature of the furnace, typically at a rate of 
3000°/sec. In fact, a gel pellet of such mass and dimensions forms a 
thick-walled sphere 356 ym diameter whose walls are opaque with carbon 
inclusions after dropping through one furnace unit (60 cm long) at 1300°C. 

A significant number of these spheres will have multiple voids, or have 
single voids with remnants of the intercellular membranes (such as shown 
in Figure 2), indicating that the glass was molten but had insufficient time 
(1 sec) to form a uniform hollow sphere (or a smooth internal surface). 

Increasing the heating time of the shell by the sequential addition of 
60 cm long furnace units at 1300° resulted in shells with decreased carbon 
inclusions and increasing diameter. With all six furnace units at 1300°, 
forming a 3.9 m long hot zone, the heating time for the shell was increased 
tc 6.1 secs; the shells formed from gel in the 3.9 m hot zone were relatively 
free of carbon inclusions. Furthermore, the average diameter of these shells 
was 407 ym, representing a 50% increase in shell volume as a consequence of 
the increased time of heating. This is surmnarized in Table 10. These 
experiments suggest that the gel-to-shell transformation is indeed rapid, 
forming a crude shell in the time (approximately 1 sec) it takes the gel 
pellet to drop approximately 0.6 m. They suggest as well that the process 
is not limited by the heat transfer rate, that the gel is indeed being heated 
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very rapidly to a fluid state. It is apparent that for gel mass as large 
as 20 ygm the rate of mass transport is relativley rapid in the gel-to-s‘'ell 
transformation, and that tia rate of glass fining is the limiting factor 
in the production of quality shells. 

Further, the rate of heat transfer by the furnace gas is sufficient to 
melt gel pat tides as massive as 200 ugm and is not the r^te-1 imiting 
mechanism of the gel-to-glass transformation. Rather, experimental 
observations indicate the gel particles as large as 200 ygm simply do not 
develop sufficient drag and buoyancy that would otherwise decrease the 
terminal velocity and thus increase the residence time of the gel particle. 
This suggests that for such large particles, the shell formation is limited 
by the large mass ana the "viscous" inertia of the gel. 

The feasibility of describing and modeling the fluid dynamics of 
the transformation, which is central to the development of a complete and 
general gel-to-shell transformation model, was examined. However, the 
complexity of the process was too formidable to model with our limited 
time and resources. Thus, rather than addressing the question of whether 
process conditions were insufficient to develop good shell uniformity, 
we chose to address the question of possible mechanisms that would affect 
the uniformity of such shells already formed. 
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5. EFFECT OF GRAVITY ON SHELL CONCENTRICITY 


5.1 Shell “Decentering ModeK A consideration in the manufacture 
of higbly-uniform glass shells is the effect of gravity on the uniformity 
of the product. The force of gravity can produce decentering of the 
interior and the exterior surfaces of a molten shell that is initially 
concentric by causing the molten glass to flow while the shell is falling. 
This motion of the interior surface can be viewed as the motion of a buoyant 
gas bubble which rises within the fluid glass shell at a rate governed 
by the buoyant force and the glass viscosity. A decentering model was 
developed which predicts the extent of this decentering of the surfaces 
of an initially-concentric shell as a function of time and the applied force. 

The motion of a small buoyant gas bubble in an infinite glass pool is 
adequately described by Stokes Law. This description is not adequate 
for molten glass snells with a thin but finite wall, wherein the "wall 
effects" must be considered, since glass flow is constrained between two 
spherical surfaces. This configuration is crudely modelled, at best, by the 
classical case of a buoyant gas bubble in a cylinder of fluid of comparable 
dimensions, such as Francis' modification of the Stokes equation [11]. Wang 
has suggested the use of the Stefan-Reynolds equation to model the effects 
of buoyant forces in molten glass shells. This treatment may be more 
satisfactory but details of the analysis have not, to our knowledge, been 
published [12]. The paucity of satisfactory models iias prompted us to 
initiate our own investigation of this question. Consequently, we have 
developed a model quantitatively describing the development of wall- 
nonuniformity in a molten shell, by gravity-driven decentering of the inner 
and outer surfaces. The model assumes that the shell initially has 
perfect wall uniformity. The decentering process is described as a function 
of time as well as glass characteristics. 

The model is based on the following conditions and assumptions: 

1) Two types of forces control the dynamic behavior of the fluid 
shel 1 : 

a) forces due to the constant acceleration field of gravity 

b) forces caused by the viscosity of the glass fluid. 

2) The fluid shell has initially two concentric, spherical surfaces 
stabilized by the surface tension of the glass 
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3) The density of the gas in the internal void is negligible ccupared 
to the glass density. 

4) The molten glass is a Newtonian fluid of constant viscosity. 

5) All stresses, other than those due to viscous forces, are 
negligible. 

6) Transient velocity distributions are negligible. 

7) The shell is thin. 

The wall nonuniformity (WNU) at time, t, is defined as: 



where 

6 is the maximum thickness of the shell at time, t, and 

TT 

is the minimum thickness of the shell at time, t. 

Since the shell is initially uniform, WNU = 0 v/hen t = 0. The mathematical 
expression of WNU developed in the model is 




where 


T _ 8paD t „ ; where p is the density of the molten 

' war 

glass, a is the acceleration, which for our case is g, the acceleration due 

to gravity;- D is the outside diamter of the shell; y is the absolute 

viscosity of the molten glass; (AR) is the aspect ratio of the shell; the 

ratio of the outside diameter of the shell to the mean wall thickness; t 

is time. The development of this model is disucssed in detail in Appendix 

IV. As an example of the model’s application consider the following 

molten glass shell and environment: D = 200 ym, (AR) = 100, y = 300 gm/cm/sec, 

2 

p = 2.4 gm/cc, and a = g = 980 cm/sec . Using equations 1 and 2 the time 
interval, t, required to develop a WNU of 0.1 is found to be 2300 sec. This 
time interval is probably 2 orders of magnitude greater than the interval 
during which the shell will be molten in its fall through a 4 m tower furnace. 
The time required to develop 1% nonuniformity (WNU = 0.01) is 260 secs. 

Thus, it is improbable that gravity effects will cause even a \% nonuniformity 
in this shell. 
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In general all three models, (KMSF, Stefan-Reynolds equation, and 
Stokes-Francis equation) indicate that gravitational decentering of shells 
is negligible for shells less than 2 mm in diameter. For our model, the 
decentering time as a function of shell diameter is shown in Figure 11; 
shell parameters are the same as above. Thus a 2 mm diameter shell will 
develop a wal 1 -nonuniformity of 0.1 in 200 seconds, which is a long time 
interval compared to a normal residence time of secs in a vertical, 4- 
meter long furnace. For a 2 mm shell, the decentering time is comparable 
to furnace residence times only for thick-walled, low aspect ratio 
shells. A 2 mm shell with 80 um walls decenters in approximately 15 secs 
to a wall -nonuniformity of 0.1. Thus gravity will significantly decenter 
only those shells having a large diameter, a long residence time in the 
furnace, a low aspect ratio, or a very low viscosity. For materials such 
as liquid metals with a viscosity near 0.01 gm/cm sec, the shell in the 
example (200 ym 00) would develop a wall -nonuniformity of 0.1 in 0.07 
seconds. The resistance of a glass shell to gravity decentering is 
essentially a consequence of its relatively large viscosity. 

We are concerned however that the model may have a serious weakness 
resulting from one of the boundary conditions employed in its development. 
Specifically, in an evaluation of the development of this model, R. S. 
Subramanian [13] has suggested replacing our assumption of no slip at the 
exterior shell surface by a condition of zero tangential stress. We agree 
with this recommendation; however, the incorporation of the zero-stress 
condition into the model will produce a much more complex model which 
probably cannot be solved in closed form. 

5.2 Preliminary Model Testing Experiments . As discussed above, the 
model for gravity-driven shell decentering predicts very low rates of 
decentering for glass shells of the small dimensions and compositions 
currently produced because the shells experience only relatively short 
residence times in the furnace. 

Testing of the model requires heating shells for much longer time 
periods than those elapsing during shell formation. One method of 
increasing heating time is to repass shells through the furnace. However, 
our past experience has shown that multiple passes of glass shells through 
a vertical furnace is rather ineffective for remelting and reforming 
the shelUeven though the sample surface is uniformly heated, the particle's 


52 



2500-1 


2000 H 


\ 

■\ 

■\ 


^ 1500- 

UJ 

</i 

LU 


1000 H 


\ 

■\ 

\ 


\ 


Figure 11. EFFECT OF SHELL DIAMETER 
ON DECENTERING TIME 


FINAL WNU = 0.1 
AR = 100 



H 1 1 1 1 r 

0 1 2 3 4 5 

SHELL DIAMETER X 10' cm 


53 


residence time and temperature history cannot be controlled easily. 

Further, air turbulence at the top of the furnace results in the loss 
of many shells to the wall of the furnace during shell injection. In 
addition, a continuous rather than cumulative re-heating of shells is 
preferred, to avoid questions concerning the effects on .i.e shells 
of repeated remelting and solidification. 

A controlled gas-jet which levitates the shell in a furnace, rarher 
than repeated passes through the furnace, is a more satisfactory method 
to increase shell residence time at its forming temperature. 

A levitator/ furnace system, with a collimated hole structure (CHS) 
levitation jet and an el liptically- focused infrared heater, was made 
available to us by Edwin Ethridge of Marshall Space Flight Center (MSFC) 
[8]. It had been designed and build by S. Dunn (Bjorksten Laboratories) 
and E. Ethridge (MSFC). We conducted a series of levitation/heating 
experiments jointly with E. Ethridge at MSFC using a batch of shells 
prepared at KMSF. The shells ranged in diameter from 500 to 1000 urn 
and were prepared from a gel of nominal composition 73% Si02» 12% K2O, 

8% Na20, and 7% 620^ (weight percent). The shells were not analyzed 
for composition, but on the basis of similar experiments, are expected 
to have considerably less alkali and boron than found in the gel. 

In general, the very large shells ('^> 1 mm) levitated easily and 
stably, requiring only a minimum of gas flow adjustments during the heating 
sequence. However, poor uniformity of sample illumination (heating), 
poor levitation stability during temperature programming and slow heating 
rates ('\>300°C/min) were significant problems encountered for heating large 
shells with this type of system. 

At 700°C, the shells generally dimpled and wrinkled in the hemisphere 
of the shell opposite the impinging gas stream, and this deformation became 
progressively worse with increasing temperature. This uneven heating 
and collapse is due to the system configuration. The shell is shielded 
from reflected radiation from below by the jet body, and partially from 
above by the thermocouple bead. Complete collapse of the shell did not 
occur. A few shells, upon reaching 900-1 000°C, turned opaque 
("frosted"), presumably due to phase separation. The shells that remained 
stable during the partial collapse did reform into spherical shells above 
1050°C, but with smaller diameters. In these shells gross changes were 
observed in the shell dimensions and surface finish. 
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The recovered shells were examined by optical and SEM microscopy at 
KHSF. Shells that were heated at 900°C or more for any length of time had 
become translucent or opaque. Optical microscopy showed a "wormy" optical 
i texture in the glass of the shells, indicative of phase separation, as 

well as large cracks and gross surface deformations. SEM microscopy 
revealed surfaces riddled with deep pitting and craters (20 to 50 pm in 
diameter), and numerous large cracks as shown in Figure 12. Whereas the 
glass surface of the unheated shells was smooth with roughness < 0.05 pm 
(except for a moderate incidence of 0.1 -0.5 pm weathering bumps), 
the remelted surface had a uniformly coarse and grainy texture (^ 0.5 pm), 
appearing almost porous, as in Figure 13. 

Shells heated to less than 900°C, while optically clear, shewed surfaces 
that were coarse and dimpled with very shallow flat-bottomed craters 
(10 to 20 pm in diameter). These craters invariably had fine diametric 
cracks, as in Figure 14. Because of the crater frequency, these shells 
resembled golf balls. The Energy Dispersive X-Ray Spectroscopy (EDXS) 
data indicated that the bottoms of the craters were depleted in alkali, with 
an alkali (Na and K) composition less than half that of the surrounding 
glass matrix. The alkali loss was localized to cratered areas; the composition 
of the surrounding matrix was similar to that of unheated glass shells, 
shown in Table 1 1 . 

The observed change in composition and degradation in shell quality, 
resulting even from mild ('v 900°C) heating, underscores the difficulty 
of remelting glass shells either to improve them [14] or to deteiTiiine 
the effects of gravity on shell WNU. These problems may be solely a function 
of glass composition; the composition may be particularly susceptible 
to phase separation and alkali loss. Other compositions may be substantially 
more stable but our experience with glass indicates a significant volatilization 
rate for alkali and boric oxide component:-; in most glasses, wht!. in the 
form of small shells. Thus the glass-melts continually change in composition 
when held at temperature for a long time, resulting in a corresponding 
increase in viscosity. Since glass viscosity is a parameter of the model 
which must be held constant (and since increasing viscosity will decrease 
the rate of gravity-induced glass flow) the model is quite difficult 
to test unless the glass composition is held constant. The change in 
composition with time may also promote phase separation, but equally 
conducive for phase separation is the relatively slow heating and cooling 
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Figure 12. Photomicrographs— Phase Separation in Shells Reheated at 1300“ C. 
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Figure 13. Change in Surface Texture Upon Reheating Sheils to 900'^ C. ! 
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Table 11. 


EDXS Elemental Analysis of Remelted Glass Shells 


EDXS Analysis 


Sample 

Treatment 

History 

Sampling 

Na/Si 

K/Si 

Control 

not heated 

averaged- 
whole surface 

0.017 + .008 

0.029 + 0.015 


1 

900°C 

whole surface 
crater bottom^ 

0.019 + .004 

0.030 + 

.016 

2 

900°C 

0.008 + .004 

0.008 + 

0.004 

3 

1300°C 

matrix 

whole surface 

0.016 + .002 
0.0092 

0.021 + 
0.001^ 

0.002 


1 Data for 5 crater and 5 matrix readings on one shell. 

2 Shell melted and reformed at 1300°C; the sole shell 
recovered from this treatment. 
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rates used in the experiments. The solution to this difficulty may be 
a) to minimize the total heating time by increasing the efficiency of 
melting the shells, e.g., with a CO 2 laser beam [15], or b) to employ 
non-volatile glass systems with very low and stable viscosity. Germania, 
a single-component glass, would be an ideal candidate on the basis of these 
criteria, except for its relatively high viscosity in the liquidus. 

O 

Although lower than the silica liquidus ('^^10 poisel, it is still higher 

than found for many alkali borosilicate and soda-lime glasses which range 
4 

from 10 to 10 poise. In testing for the effect of gravity on molten shells, 
such viscosities would require unacceptably long sample residence times in 
a vertical furnace, or higher temperatures than are currently possible 
with the CHS levitator and radiant furnace described above. 

Since the viscosity of pure germania glass is considerably higher than 
desired for reheating experiments, binary glasses in the silica and 
germania glass systems were investigated. Phase diagram data eliminated 
the alkaline earth silicates as candidates. Even though this system has 
the advantage of a stable composition due to the non-volatility of the 
alkaline earth oxide glass modifier, it has the problems of a) phase 
separation in the liquid state due to immiscibility over a large concentration 
range, and b) an excessively high eutectic temperature (> 1400°) [16]. 

Alkali silicates were undesirable since high concent ations of alkali 
('\^15-25 mol percent) are required to significantly reduce the liquidus 
temperature (to 1000-1200°C) [16]. Such high concentrations of alkali 
render the glass susceptible to alkali volatilization and the consequent 
composition and viscosity changes upon reheating. 

The binary germania system has an inherent advantage. The liquidus 
temperature of pure germania is already in the desired range. The aodition 
of moderate amounts of alkali oxides, e.g. mole percent, will bring 
it well within the desired range, from 1150°C to 900-100u°C; the eutectics 
occur at 4-10 mol percent M 2 O [17]. The alkali ai^e desired solely for 
viscosity modification, and in the germania system, small amounts of alkali 
oxide, such as Na20, very sharply reduce the viscosity. For example, 
the addition of 1.2 mole percent Na 20 to pure Ge 20 reduces it to 20 poise 
[17,18]. The addition of lithium and potaasium oxides yield similar 
results. While high alkali oxide ('^5-7 mol percent) germania glasses are 
expected to be moderately volatile and compositional ly varying at the 
liquidus temperature, the correspondingly low glass viscosity should require 
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a very short sample residence time in the furnace (during experiments 
testing for gravity effects on the molten shell). Because of the low 
operating temperatures and the initially low concentration of alkali, low- 
alkali oxide germania glasses (< 1 percent) are expected to have only 
negligible alkali loss rates during the course of the reheating experiments, 
and the composition and viscosity should remain constant in time. Thus 
a series of glasses, with a wide range in viscosities, may be possible for 
testing the gravity-driven decentering model. 

Although of poor quality (wall uniformity) germania shells were made 
by the metal-organic gel method several years ago at KMSF using precipitated 
powder from the hydrolysis of germanium ethoxide. More recently, the 
synthesis of homogeneous germania gels was attempted by hydrolysis of ethanol ic 
solutions of a) germanium tetrachloride, b) germanium tetrachloride acidified 
with oxalic acid [19], and c) germanium tetraethoxide [20,21]. However, 
these reactions yielded either crystalline or gelatinous precipitates 
which dried to a fine, chalky powder rather than a glassy friable gel. 

These gels were also devoid of residual organics, as confirmed by TGA 
and PGC. 

Superior homogeneous gels of germania and the more reactive sodium 

germanate (5 mol % Na20) were synthesized by controlled partial hydrolysis 

of alcohol solutions of germanium and germanium/sodium alkoxide mixture 

(80% of stoichiometric amounts of H2O) at -50°C. These reactions yielded 

homogeneous, translucent gels which, when dry, were glassy and friable. 

Sodium german* te shells as large as 600 ym have been produced at 1300°C 

from the freshly-prepared corresponding gel. These gels are very sensitive 

to atmospheric moisture. They hydrolyze readily and completely, with 

consequent loss of all the residual organics and a conconitant reduction 

2 2 

of the high surface area (from 460 m /gm to less than 10 m /gm). Unlike 
the silicate gels, this hydrolyzed sodium genmanate gel produced only small 
thick-walled shells of very poor uniformity. 

A possible approach to get around this problem is the incorporation 
intothe gel of an auxiliary blowing agent which is inmunc to hydrolysis. 

Urea, a compound used in blowing of silicate based glass shells, was 
successfully incorporated into the germania gel in an apparently homogeneous 
manner, but shell blowing experiments could not be done before the expiration 
of this contract. 
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6. CONCLUSIONS 

The purpose of this work has been to study and model the mechanisms 
that produce glass shells from metal -organic gels, and thereby to provide 
a basis for an assessment of the benefits of manufacturing shells in space. 
The objectives have been to determine: 1) the eftccts f gel pyrolysis, 
and thus the nature of the gel-to-shell transformation, 2) the effective 
heating rate of a gel particle falling through a vertical furnace, and 
3) the effect of gravity on the shell concentricity (i.e., shell-wall 
uniformity) as the molten shell drops through the furnace. 

A gel particle falling through a vertical furnace is rapidly pyrolyzed 
and transformed to a hollow shell. For the size of gel currently used 
('V 20 ygm), that transformation occurs within one second (@ 1300°C) and 
is rapid with respect to the time needed for complete fining of the 
glass in the shell. The transformation to a shell is comprised of several 
distinct processes: a) surface closure of the porous structure of the 
gel, either by low-temperature (< 400®C) sintering or else by hydrolysis, 
b) formation of a closed-cell foam structure ('\^ 400°C), c) expansion and 
spheridizatior of the particle ('\- 600-800°C), d) coalescence )f the foam 
to form a single-void shell ('v- 800-1000°C), and e) fining ot che glass 
walls of the shell. These processes have been observed during gel 
pyrolysis by hot stage microscopy, or have been inferred from pyrolyzed- 
gel intermediates recovered from shell-forming experiments. The driving 
force for the formation of the initial closed-cell foam structure in the 
gel is supplied by the gaseous pyrolysates generated from the gel. These 
gases are typically comprised of H 2 O, CO 2 , ethylene, ethane, and methane, 
and are generated from the metal alkoxides, the si lands and the bicarbonates 
in the gel at fairly low temperatures. The pyrolysis decomposition of 
the metal -alkoxides also produces elemental carbon, resulting in a 
characteristic black appearance to the gel. The oxidation of this carbon, 
as well as the decomposition of metal carbonates in the gel, produces the 
CO 2 that drives the further expansion of the foamed spheroid and the 
final shell. Undoubtedly, most of the pyrolysates are rapidly lost and 
may be ineffective in shell -blowing, except to separate a charge of gel 
powder into a w«»ll -dispersed plume upon injection into the furnace. 

The average metal -organic gel loses more than 20% of its mass as gas, 
while much less than 1% rema'-ns to form the final shell. These two 
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conclusions are supported by two observations. Gels that were pyrolyzed 
to remove all of the gaseous pyrolysates ('\/ 26 % of the initial gel mass) 
formed shells in yield and quality similar to those Tiade from unpyrolyzed 
gel. In addition, metal-organic gels that were hydrolyzed (under argo 
or in air) to remove all of the organic species (except for tr'ice levels) 
produced shells of similar yield and quality as those from unhyot-olyzed 
gels. In the latter case, the blowing qas, C02» is derived primarily from 
carbonates in the gel, and comprises 1% of the total gel mass; yet only 
.01% ends up among the final blowing gases in the shell. (The other 
constituents of the residual gases in the shells are Og, N 2 and sometimes 
H 2 O. The O 2 and N 2 permeates the shell from the furnace ambience.) 

The various .asorption and decomposition reactions that occur during 
gel pyrolysis constitute a thermal inertia, or latent heat, in the gel, 
and increase the specific heat of the gel over that of the corresponding 
glass. However, the oxidation of the organic pyrolysates and the reduction 
in surface free energy (due to sintering) are exothermic and partially 
compensate for the endothermic reactions. ;'he average specific heat for 
unhydrolyzed metal-organic gel is less than twice ('\^ 1.5-2 j/gm K) that 
for glass. In any event, this latent heat has no effect on shell production, 
as was evidenced by the shells produced from pyrolyzed gels. 

The results of the heat-transfer modeling and the experiments in 
support of this modeling indicate that gel particles as large as 20 pgm, 
which yield 400-600 ym shells (depending on furnace conditions), are 
heated very rapidly, at a rate of approximately 3000°C/sec. The gel- 
to-glass shell transformation is completed in 1 sec. The final fining 
of the glass, to reinove the remaining carbon ano gas inclusions in the 
wall, takes considerably longer depending on tlie furnace gas and 
temperature (e.g., >5 secs at 1300°C, in air). Gel particles j. large as 
200 ygm are also heated i-apidly to the temperature of the furnace ambience, 
at o rate of approximately 2000°/sec. However these massive particles 
develop a high terminal velocity and fall deeply into the furnace (> 1 m) 
beforv’ reaching that temperature. Because of the larger mass, the gel 
transformation is significantly slower, and the particles do not expand 
rapidly enough to generate sufficient buoyancy and drag to decrease the 
terminal velocity 0 ^ the particles. Consequently such massive gel panicles 
have a very short residence time (< 2 secs) and cannot transform completely 
into shells. The crux of this problem is the high terminal velocity. 
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of the particles. Thus the problem can be solved only by increasing 
the furnace length, by reducing the terminal velocity with counter- 
flowing gas (or gas-jet levitation); or by reducing the gravitational 
acceleration on the particle. 

The model for gravity-induced decentering of shells predicts 
significant decentering rates for moUen glass shells less than 2 mm 
in diameter. There is a real concern that the model has improper 
boundary conditions for the external shell surface. With the apparently 
appropriate boundary conditions extremely rapid and extensive decentering 
of the internal and external shell surfaces are suggested. However, 
such extensive decentering has not been observed in current shell 
manufacture. These observations suggest that either these boundary conditions 
are not correct, or an active concentering mechanism exists to counteract 
the effect of gravity. The latter appears to be the more likely answer. 

Efforts to test the validity of the decentering model have been 
hampered by problems related to uniform and controlled shell heating, 
stable levitation, and stability of the glass; compositions. Consequently 
these empirical studies have not yet been completed. 
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7. RECOMMENDATIONS 


We believe that a significant advance in the understanding of the 
processes involved in glass shell synthesis from metal -organic gels 
has been achieved from work done under this contract. On the basis of 
this knowledge, plus the additional understanding which has accrued from 
related work under our Department of Energy contract over the past three 
years, we can identify incomplete and unstudied areas for further research. 

We have been motivated in our research by the ultimate goals of 1) understanding 
the mechansims which control the formation of highly uniform glass shells 
suitable for inertial fusion energy research and 2) preparing for the 
potential necessity of making in space large shells which may not be 
obtainable in a terrestrial laboratory. Theses goals also motivate our 
suggestions and recommendations for fu'^ure work. These recommendations 
can perhaps best be made in the context of further mechanistic studies 
required for three processes that affect the formation of uniform glass 
shells from gel. Each part of any future research should contain both a 
mathematical and an experimental component. 

7. 1 Mass Transport Study 

Mathematical Model. Develop a second generation heat transfer/particle 
trajectory/ fluid mechanical model for gel particles forming a shell in 
a vertical furnace. This model should be able to handle inhomogeneities in 
physical and chemical properties for the gel and shell during the shell 
forming process. 

Experimental Program. Conduct experiments to collect thermodynamic 
and chemical kinetic data required for mathematical modeling. 

A successful model would predict the time and temperature requireo to 
form a complete bubble-free shell from a given size gel particle of known 
composition. 

7.2 Shell Concentricity Study 

MathemcMcal Model. Experimental evidence from shel 1 -forming and 
gel oyrolysis studies suggest the existence of an active concentering 
mechanism(s) , that produces shells with uniform walls from asymmetric spheroids. 
Among the concepts that have been suggested and/or partly investigated we 
propose that the following concentering forces be studied in more detail, 
a) .he tumbling action of the shell, b) surface tension gradients in the shell 
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surface caused by compositional gradients resulting from differential 
vaporization of alkali and boron oxides, c) viscosity gradients resulting from 
the composition (or temperature) gradients, or d) expansion-contraction 
cycles caused by applied force or by heating and cooling cycles due to small, 
localized temperature gradients in the furnace. 

Experimental Model. Conduct experiments to test the validity of models 
as they are developed. The role of composition gradients (b and c above) 
should be examined by comparing shell blowing experiments using gel 
starting materials with a wide range of volatile components and variations 
in chemical homogeneity. 

7.3 Shell Decentering Study 

Mathematical Model. The current gravity- induced decentering model 
(of an initially concentric shell) may have serious weaknesses due to the 
boundary conditions employed. We propose the development of a model based 
rigorously on the principles of Low Reynold's Number Hydrodynamics with 
appropriate use of computers to achieve numerical soluti ns. 

Experimental Model. Considerable effort has been expended in preliminary 
experiments concerning the effect of gravity on the concentricity of the 
inner and outer shell surfaces. The progress of this work has been hampered 
by problems of a) compositional instability of the molten glass shells and 
b) inadequate methods for containerless reheating of characterized shells. 

These problems completely mask the more subtle effects of gravity. We 
recommend the continuation of this line of research, with emphasis on the 
development and use of a) glass shells having a suitably low viscosity, 
a stable composition, and resistance to p>.ase separation and devitrification, 
b) a containerless method for uniform, controlled heating of the shell surface, 
such as an acoustic or air-jet levitation with a C 02 ~laser heating system, 
and c) shell -forming and shell -reheating experiments aboard KC-135 flights. 

In this manner, the uniformity of the shell walls could be investigated 
as a function of gravity and heating time. 
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APPENDIX I 


Heat Transfer/Thermodynamic/Ruid Mechanical Model of Processes Occurring 
During the Manufacture of Glass Shells in Vertical Furnaces^ 

The purpose of this appendix is to outline the development of a first- 
generation model of the processes occurring during the manufacture of 
glass shells in vertical furnaces. The purposes of the model are: 

1. To investigate the heat transfer and fluid mechanical 
phenomena associated with the manufacture of glass shells 
in vertical furnaces. 

2. To delineate the furnace parameters that are important to the 
manufacture of high-quality glass shells, and their quantitative 
influence on the quality of the product. 

3. To predict the largest shells of acceptable quality that can 
be manufactured in a 1 g environment with existing equipment. 

4. To specify furnace equipment that will be required to manufacture 
glass shells for future requirements, which will have more 
stringent specifications. 

This type of model is primarily the embodiment of some of the fundamental 
principles of classical physics. Effective use of the model requires the 
availability of additional information, such as materials property data, 
and insight to the processes resulting from specific empirical investigations. 

This first generation model was developed at the beginning of the 
program when little information was available on the details of the processes 
occurring in the shells and their precursors. One consequence is the 
assunption that, at every instant, the shells and their precursors are 
homogeneous. During the course of the program signivlcant insight has 
been gained regarding the temperature ranges in which specific physical 
and chemical processes occur, and relevant quantitative materials properties. 
Consequently, a second generation model is currently being developed in 
which the shells and their precursors are not assumed to be homogeneous. 

The furnaces with which we are concerned have interior heated regions 
whose walls are vertical circular cylinders. The cylinder diameters are 
about 3 inches, and the lengths of the heated sections have varied from 

♦Appendix I was authored by R. 3. Jacobs of R. B, Jacobs Associates, 
consultant to current contract. 
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2 to 12 feet; about 2 feet of unheated region ace at the bottom end. The 
furnaces are heated by electrical resistance heaters, usually wired so 
that several axial sections can be heated and their temperature controlled 
independently. In existing equipment, the operating temperatures are 
varied from 800®C to 1500® C. The heated region is usually filled with 
static, atmospheric air; the only flows are those due to natural convection 
and those that may be caused by the presence of material associated with the 
manufacture of the glass shells. The species of gas within the furnaces 
can be changed, and generally axial flows can be imposed upon the gas. 

The raw materials from which the glass shells are made are classified 
as gel powders and pellets. (For convenience in the following discussions, 
the raw material will be called "frit"; it need not be frit, but may just 
as well be cast, extruded, etc.) The fnt is introduced, at room temperature, 
into the upper end of the furnace, and falls vertically downward (except 
for possible effects of radial and tangential components of the gas velocity, 
and of the method for introducing frit). As it falls, the frit and its 
descendents are heated by the furnace walls and the furnace gas. Heat 
transfer involves both radiation and convection, the latter including gas 
conduction. Additionally, the frit and its descendents interact mechanically 
with the furnace gas; this interaction, along with the effects of gravity, 
determines the motion of the particles. As a result of the heating, the 
constituents of the frit begin to react to form glass and a blowing gas. 

Under proper conditions, the physical and diemical processes produce high- 
quality, thin-walled glass shells. The products of these processes are 
collected at the lower end of the furnace. 

It is postulated here, and can be inferred from experience, that the 
manufacture of high-quality glass shells requires that the gel material 
be heated to a sufficiently high tempe'^ature and maintained there for a 
sufficiently long time interval. We may therefore conclude that the rates 
of heat transfer from the furnace to the materials therein, and the residence 
times of the materials within the furnace are of primary importance. 

Therefore, our primary objectives here are to investigate the phenomena and 
variables that influence the heat transfer and residence times, and to 
relate these to characteristics and properties of the gel and glass shells, 
and to design and operating parameters of the furnaces. 
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The remainder of this development therefore consists of two 
sections: 

1. Heat transfer considerations which embody the requirements of 
the First Law of Thermodynamics. 

2. Particle motion considerations which embody the *^quirements of 
Newton's Second Law. 

The model is expressed here in differential form and has also been 
expressed in finite difference form for programing for parametric studies. 

HEAT TRANSFER 

The heat transfer from the furnace (walls and gas) to a particle 
therein (gel particle or shell) is described by: 

f =h(Tg-T).aE(Tg‘'-T''). 0) 

The left-hand term is the rate of increase of energy of the particle 
being heated. The middle term describes the rate of convection heat transfer 
(including gas conduction) from the furnace gas to the particle. The 
right-hand term describes the rate of radiant heat transfer from the furnace 
walls and gas to the particle; we assume that the walls and gas are 
effectively at the same temperature. 

Equation i is applied to individual particles from the time that they 
are introduced as frit to the time that they exist as shells. We ignore 
mechanical and thermal interactions between the particles themselves, and 
mechanical interaction between the particles and the walls. It has also 
been assumed, in the derivation of equation 1, that the temperature of the 
particle is always uniform. This assumption, together with the assumption 
of material homogeneity, give rise to the most severe limitations of this 
first generation model. It is justified, from a thermal history viewpoint, 
by the fact that the time intervals required to equilibrate temperature 
differences within the particles are short relative to the time intervals 
associated with the heating of the particles. The elimination of the 
temperature-uniformity and material homogeneity assimiption is the basis 
of the second-generation model mentioned above. 

The variables and parameters in equation 1 are: 
m is the mass of the particle. As the manufacturing process proceeds, 

m changes. If these changes v/ere known as a function of temperature. 


72 



they could be readily managed in the development of solutions to 
equation 1; however, this was not the case. We did know that a 
significant mass (of the order of 15%) was lost below 600°C as the 
result of the loss of blowing gas. Therefore, in this model it is 
assumed that there is a single, instantaneous loss of mass. This loss 
occurs when the particle attains the temperature T-j ; the fraction of 
the initial mass which is lost is (1 - r). T^ and r can be treated as 
parameters whose values are determined by efforts to obtain agreement 
between the model and empirical information. On the basis of early 
experiments, initial values for these parameters can be T^ = 500°C and 
r = 0.85. 

C is the effective specific heat of the particle being heated. It should 
account for the various latent heats associated with the manufacturing 
process; e.g., heats of chemical reaction and phase transitions. There 
was also a dearth of information in this area. 

Ag is the surface area, of the particle, that is exposed to the heating 
environment. It varies with time, and may be correlatable with 
temperature. However, it is not now possible to express A as a 
function of time or temperature. In this model, we specify a geometry 
for the initial frit particle: a right circular cylinder of diameter 
and length (2 D^). This limitation can be eliminated as discussed 
below. We assume thatat aspecified temperature (T 2 ), the frit particle 
instantaneously transforms to a sphere of diameter, D, the final diameter 
of the shell. As will be discussed below, and D will be related by 
a parameter which, along with T 2 might be determined by correlating 
empirical information to fit the model. 

T is the absolute temperature of the particle being heated. 

t is time. 

h is a heat transfer coefficient that characterizes the convective (and 
gas conductive) heat transfer from the furnace gas to the particle being 
heated. It depends upon the geometry of the particle, properties of che 
furnace gas, and the flow field of the furnace gas relative to the 
particle. For this model, we express it in terms of fairly reliable 
steady-state theoretical formulations and empirical correlations. This 
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is partially justified because the time intervals required to develop 
steady-state flow fields are short relative to the residence times 
in the furnaces. The major uncertainties stem from lack of information 
regarding the geometries of the particles as they pass through the 
furnace. The same assumptions are used to assist in the calculation 
of h as were discussed in connection with A^. 

As we shall see later, in our discussion of residence times of 
particles in the furnace, the Reynold's numbers associated with the 
particles in the furnace are small. From McAdams^ we find that we 
shall be somewhat conservative (i.e., predicted temperature changes 
are too small) if we use: 

k 

h = 2 where (2) 

kg is the thermal conductivity of the furnace gas, and 

L is a characteristic dimension of the particle. 

For the frit particle, L is assumed to be the aritnmetic average of the 

length and diameter; for the shell, L is its diameter, D. 

Tg is the absolute temperature of the furnace gas (and walls). In 

general, Tg is a function of position in the furnace, and therefore, 
in equation 1, a function of t. In our experiments, this function is 
known and is therefore used in developing solutions to our model. 

-12 

a is the Stefan-Boltzmann constant, and is approximately 5.67 x 10 
watt/cm^K^. 

E is the effective absorptance of the particle to the radiation emitted 
by the furnace walls and gas. It is a function of the geometries 
of the furnace and particle, and of properties of the furnace walls, 
furnace gas, and particle. Because the object being heated is tiny 
relative to the size of the furnace, E is effectively equal to the 
absorptivity of the particle to the radiation emitted by the furnace. 

It is probable that the frit absorbs essentially all of the radiation 
incident upon it, and that the clear glassy material will absorb a 
negligible amount of the radiation incident upon it. We assume, 
therefore, in this model that E = 1 until the temperature of the particle 
becomes T^; at this time the value of E instantaneously becomes zero 
(E = 0). As indicated with T^ and T 2 , T^ can also be used as a 
correlating number. 
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We are now in a position to formulate the preceding discussion into 
relations that can be substituted into equation 1. The results are: 


For T < , 


27^ r 





(Tg - T) * - t'*). 


For < T 4 : T 2 , 


‘’ql '^f 
2.5 



n 


dt 



D 


J/3 


(Tg - T) 


+ a(T 



(3a) 


(3b) 


For* Tp ^ T ^ 3 * 


Pgl ^ ^ ‘^9- 


T^, 


For T > T^, 

P,1 f 2 ^ - T). 


(3c) 


(3d) 


Equations 3a, 3b, 3c, and 3d (i.e., equations 3) are the final embodiments, 
in differential form, of our energy equation for the temperature ranges 
of interest. The new symbols which have not yet been defined are: 

= is the mass density of tne glass in the shell. 

* is the effective specific heat of the frit material. 



= is the mass density of the frit material, 
n = is the thickness of the shell wall. 

C 1 = is the effective specific heat of the glass in the shell wall. 

9' 

It is reasonable to assume that, for a specific shell manufacturing 
process, p^, p^-j , and r should be fixed. Thus, from the above definition 
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of A, there should be a unique linear relation between (D/D^) and (D/n). 

In order to become numerically oriented, we have used some of our data to 

\l/3 


calculate 



We obtained that 


_D_ 

D, 


» .16 


There is no reason we believe that this relation is valid for any situation 
except for the one from which it was derived. 


PARTICLE POSITION 


The purpose of this development is to derive relations that will allow 
us to predict the position of a particles, as a function of time, as it falls 
through the furnace. By using these relations, in conjunction with equations 
3, we can predict the temperature of the particle as a f unci ton of its 
position in the furnace. 

The position of the particle, as a function of time, is computed by 
integrating its velocity with respect to time. The velocity of the particle 
is obtained by applying Newton's Second Law to the particle falling through 
the furnace gas under the action of gravity. The effect of the gas upon the 
velocity of the particle is manifested through the "drag" force exerted 
by the gas on the particle. As this depends upon the geometry of the particle, 
and as the mass of the particle enters directly into Newton's Second Law, 
the computation of particle position is divided into four temperature regimes, 
as was the computation of the particle temperature. 

There are three forces acting on the particle as it falls through the 
furnace: drag, buoyancy, and gravity. The drag force may be represented 

satisfactorily by Stokes' Law. The resultant differential equations for 
the velocity (U) and the displacement (S) of the particle are: 


and 


f = (y*6U). 

(4) 

f 

(5) 

/ V P \ 


^ “ 9 6 - -15^ . 

(6) 
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and 


12.5 u A 

6 . (7) 

D m 

g = acceler'* cion of gravity. 

V - volume of particle. 

Pg = density of the furnace gas. 

m * mass of particle. 

y- = absolute viscosity of furnace gas. 

A = frontal area of particle. 

IT = characteristic dimension in Reynolds' :juinber. 

The deta'led formulations of the parameters in equations 6 and 7 
depend upon the temperature ranges discussed previously. For each temperature, 

Y and 6 can be taken as conscant {if average values be used for and Pg). 
Then, equations 4 and 5 can be easily solved for U and S: 


(y + 6U) = (y + 6Uq) e * ^0^ (8) 



U = Ug and S = Sg when t * tg. 


Because of the non-linarity of equations 3, and because of the variable 
properties therein, equations " are solved numerically by means of computers. 

H is therefore convenient to solve the particle position relations lumerically 
also. As a result, our working model uses equations 4 aid 5, rather than 
equations 8 and 9. 

Ii. producing our previously discussed concepts into equations 6 and 7, 
we develop our «>xpressions for y and 6 for ou- various temperature regi.nes. 



77 



5 



For T-j < T 4 T 2 



For T > (contains < T < and T > T^), 



AR (the aspect ratio) is the ratio of the diameter of the shell to its 
average wall th-'-kness. It should be noted that the foregoing model is 
expressed in tL of furnace characteristics, properties of the materials 
involved in the manufacturing process, and geometrical characteristics 
of the shell. Summarizing, our model is described by; 

1. Equations 3a, 4, 5, 6a and 7a prior to the mass loss, and 
characterized by the parameter (1 - r). 

2. After the mass loss, but prior to the transformation from rrit to 
shell, equations 3b, 4, 5, 6b and 7b. 

3. After the transformation from frit to shell, but prior to becoming 
transparent to radiation, equations 3c, 4, 5, 6c and 7c. 

4. After becoming transparent to radiation, equations 3d, 4, 5, 6c, and 
7c. 

It should be noteo, that the foregoing equations were developed 
so that the geometry of only the shell appears explicitly; the geometry 
of the gel particle does not appear. This was accomplished by assuming 
t.idt all gel particles are geometrically similar; namely, the length 
of the particles is twice their diameters. This is, of course, usually 
not the case The equations can be modified in the following ways to 
eliminate this limitation so that the r^te of heating is explicitly a 
function oi no pellet dimensions of diameter (D) and length (L), each 
of which in this case is an independent parameter. 

Thus, equation (3a) is changed to: 
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51 

dt 




(“fi 


(Tg - T) ^ o(Tg^ 



Equation (3b) is changed to: 


(3ab 


D^2 ^fZ f^f2 ^f2 £T 

+ '-fo) 

2 


^ *^g 

( 0^2 ^ ^fZ ^ 


(Tg - T) + o(Tg'> - t'*). 


Equation (3c) and (3d) are unchanged. 

Subscript f refers to frit (or gel) particle. 
Subscript 1 refers to particles when T .$ . 

Subscript 2 refers to particle when < T ^ T 2 . 
D is particle diameter. 

L is particle length. 

Equation (6a) is changed to: 

Y = g( 1 - 

Equation (6b) is changed to: 

Y = g(i - -^). 

Pf2 

Equation (6c) is unchanged. 

Equation (7a) becomes 

^ ^ _ 10_ ^ f °fl ^1 . 

9 (Dfi + L^^) 

Equation (7b) becomes 

^ ^ M ^ 4 ^f2 *~f2 . 

^ (Df2 + '-f2) ^fz ^fZ Pf2 


(3b!) 


(6a^ 


(6b^) 


(7a^ 


(7b^ 
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Equation (7c) is unchanged. 

Ug is absolute viscosity of furnace gas. 

P is mass density. 

References : 

1. McAdams, W. H. Heat Transmission. New York: McGraw Hill, 1954. 
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APPENDIX II* 

GELSHEL Program 

A program (GELSHEL) has been developed to calculate the position, 
thermal characteristics and the temperature history of a particle dropped 
through a vertical furnace. The governing differential equations are 
discussed in Appendix I and will not be repeated here. The program utilizes 
the forward finite difference of approximation to the differential equations. 

The program keeps track of the particle as it undergoes four temperature 
dependent transformations shown in Figure 1. The particle is dropped from 
rest into the top of the furnace. Initially, the input particle is at 

room temperature,, its mass and IR absorptance are m^ and a-j , its geometry 

is that of a right circular cylinder of diameter and length , L 2 and its 
specific heat is (1 + X^)Cp The additional X is used to approximate 
endothermic reactions that take place in the particle. 

Upon reaching T^ the particle undergoes changes in mass (Mj ‘I 2 ), 
dimensions (D^ -► D 2 , L^ specific heat (X^ -*■ X 2 ). At T^, the 

endothermic reactions are deemed complete and the particle's specific 
heat reverts to Cp,the specific heat of glass at T 2 , the geometry changes 
from that of a right circular cylinder to that of a spherical shell. The 
final transformation occurs when the particle temperature reaches T^. At 
this temperature, the IR absorptance changes from to The execution 

of the program is terminated when the particle exits the furnace (S = 370 cm) 

The program accepts the following input data: 

M 2 and of particle 

^^2’ *^3’ ^1 ^2 pa»*ticle 

f^glass’ h* h' °‘2 
^1 » ^2* T^, T^ 

Furnace Temperature Profile (T^ vs. S) 

Furnace Gas Code (1 ,2, 3, 4, 5, 6, 7) 

Typical values for an input data set might be as follows: 

ly M^(g): 19.2E-6, 13.8E-6, 13.8E-6 

D^, D 2 , D2(cm): 266E-4, 333E-4, 396E-4 

L^, L2(cm): 2b2E-4, 365E-4 

♦Appendix II was authored by R. J. Siinns of KliS Fusion, Inc. 
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PgUss<3''“> = 

2.4 


CM 

X 

X 

o 

« 

o 

o 

• 

< 

0^1 ,02 • 

1.0, 0.0 


Ti.T2.T3.T4(C): 

500, 1000, 850, 1250 


S(cm): 

0, 4, 7, 10.5, 18, 21 


T^(C): 

750, 1000, 1200, 1300, 1400, 1425 


S(cm) : 

25, 31.5, 46.5, 54, 370 


Tf(C): 

1450, 1475, 1500, 1510, 1510 


Furnace Gas : 

1 (air) 


The specific heat (Cp) of the particle is assumed to be that of glass 
which is temperature dependent. A table of the specifi: heat of glass 
vs. temperature is contained in arrays within a subroutine. Linear 
interpolation between table values is used to update the specific heat 
of the particle at each interaction. 

Also contained in arrays are the pertinent temperature-depencent 
properties (viscosity, conductivity and density) of seven furnace gases: 
air, argon, carbon dioxide, nitrogen, oxygen, steam and helium. Once the 
furnace gas is selected, the program writes a working table of the temperature 
dependent properties. Again, linear interpolation is used to update the 
viscosity (y^), thermal conductivity (k^) and density (p^) of the furnace 
gas at each interaction. 

At the start, the temperature of the particle (Tp) is 20°Ci,and its 
d'’>placement (S) from the top of the furnace and its velocity are zero. 

Using the value of S, the program determines the temperature (T^) of the 
furnace gas from the temperature profile table with the value of T^, the 
program then determines the furnace gas viscosity (y^), thermal conductivity 
(k^) and density (p^). Using Tp, k^, and T^ the program calculates the . 
temperature rise of the particle from its thermal capacitance, the rates 
of heat transfer to the particle by convection and radiation and the time- 
step. It then uses the nia:s and dimensions of the particle and the values 
of y^ and p^ to calculate the "new" displacement of the particle. This 
process is repeated until the particle exits the furnace. 

The output of the program, in addition to the input data, consists 
of the following: 


V A 
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1. Working tables of temperature-dependent properties. 

2. Temperatures, times and displacements at which the particle 
transformations occurred. 

3. Particle temperature, displacement, velocity and specific 
heat histories. 

4. Values of selected expressions used in the calculations of (3) 
above. 

5. Plot of the particle temperature history. 
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Fig-ire 1. Temperature-Dependent Transformation of 

Gel-Particles Dropping Through Glass furnace. 
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APPENDIX III* 




Radial Temperature Difference Model 


A program, TGRAD, has been written to calculate the radial temperature 
profile history of a spherical gel particle. The gel particle was divided 
into eight equal -mass sub-volumes, or nodes. The inner node is a sphere 
and the remaining seven nodes are spherical shells. For nodes internal 
to the gel particle the program utilizes the forward finite difference 
approximation to the one-dimensional differential equation of heat 
conduction. 




r T 1 - T T - T 1 

n-1 n ^ n+1 n I 

. Vl, n n + 1 j 


where 

T^' = Temperature of node n at the end of the time-step. 

T^ = Temperature of node n at the beginning of the time-step. 

At = Time-step. 

= Thermal capacitance of node n. 

R ^ ^ = Thenmal resistance between nodes n-1 and n. 
n-1 ,n 

The thermal capacitance of a node is given by 

C = mCp(T) 

where 

m = Mass of the node. 

Cp(T) = Specific heat of the nodal material at the temperature T. 

The thermal distance between adjacent nodes in the form of spherical shells 
is given by 



2it Dg' D.' k(T) 


♦Appendix III was authored by R. J. Simms uf KMS Fusion, Inc. 
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where 


D' = "Average" diameter of a shell, i.e., divides the shell into i»>u 
shells of equal mass. 

Dq* * "Average" diameter of outer-shell 

D^. ' = "Average" diameter of inner shell 

k(T) = Thermal conductivity of the nodal material at the temperature 
T, w/cm *°C. 

The boundary condition' is that the heat input to the outermost shell is by 
convection and radiation; this is calculated from the same heat transfer 
relationship used in the GELSHEL code (described in Appendix II). Initially, 
the gel particle is at room temperature. 

The specific heat and thermal conductivity of the gel particle are 
parameterized by linear approximations of the type A + B * T to published 
data. For specific heat, glass properties are used. For thermal conductivity, 
two extremes are used: glass properties for tightly compacted material 

and air for loosely compacted material. 

The temperature of the furnace and furnace gas is a function of the 
distance from the top of the furnace. Hence, as the gel particle descends 
through the furnace, the ambient temperature changes. The rate of descent 
of the gel particle then determines the apparent rate of change of the furnace 
temperature. For short time periods of the order of 20 msec, a 20 pg 
gel particle at the top of the furnace descends only about 2 millimeters 
which causes an ambient temperature change of only 12°C. Hence, the furnace 
and gas (air) temperature was held constant for any given problem and the 
program was executed for just 20 to 50 msec problem time which was long enough 
to determine the maximum radial temperature profile of the gel particle. 

For a furnace temperature of 1500°C, the maximum radial temperature 
difference in a 20 pgm, 296 ym diameter gel particle is 80°C for 
tightly compacted material and 1390°C for loosely compacted material as 
shown in Figure 1. For a furnace temperature of 750°C, the maximum 
temperature difference is 22°Cfor tightly compacted material and 674°C for 
loosely compacted material. For a 180 pgm, 608 pm gel particle, the maximum 
radial temperature difference is 105®C, occurring 15 msec after injection 
into a 1500°C environment (See Figure 2). 
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If the furnace environment is assumed to be 750'®C, (shown in Figures 
3 & 4} a more realistic condition, a maximum radial temperature difference 
of 25®C develops in the same particle 22.5 .nsec after injection. In all 
of these cases, the gradients maximize within approximately 20 msec and 
diminish rapidly with time. Consequently, one can safely make the 
simplifying assumption of a negligible temperature gradient for such gel 
particles in heat transfer calculations. 
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APPENDIX IV 


A Model for the Prediction of Wall Non-Uniformity, in Thin, Liquid. 

Spherical Shells, Resulting from Acceleration (e.g.. Gravity) Fields .* 

The purpcse of this analysis is to develop a model that gives insight 
to the wall non-uniformity that is produced in a thin, liquid, spnerical shell 
by gravity. A primary guideline for this development is that the analyses 
be extremely straight-forward; although it is feasible to develop very 
accurate models for this type of problem, such a project would be complex, 
long-range, and expensive. 

The general notation is illustrated in Figure 1. Although it need not 
generally be the case, we assume that the liquid shell is bounded, both 
internally and externally, by spherical surfaces. This assumption is 
justified because we are primarily interested in shells of molten glass, 
for which the effects of surface tension are relatively large. The spherical 
region inside of the shell is occupied by a gas with a density that is 
negligible relative to the density of the liquid comprising the shell. 
Additionally, the region exterior to the shell is occupied by a gas: e.g., 
a furnace gas. Because of the very large viscosities of the molten glasses 
in wfiich we are interested, we assume that two types of forces control 
the dynamic behavior of the shell: 

1. Body forces due to a constant acceleration field, e.g., gravity. 

2. Forces caused by the viscosity of the liquids. 

These assumptions place our model into a rather well-defined, very-important 
region of fluid mechanics, known as Low Reynolds Number Hydrodynamics or 
Creeping Flow. Although these types of problems have been studied intensely 
for over a hundred years, the problem in which we are interested has not 
been solved, and, as indicated above, cannot be solved within the guideline 
we have been working. 

We therefore have adopted the following simplistic image. Consider that, 
as some time t (say, t = 0), the interior and exterior surfaces of the shell 
are concentric. Immediately, the acceleration field causes the gas bubbles 
to "rise" within the liquid shell; or, conversely we may visualize that the 
liquid flows "down" the narrow space between the interior and exterior 
surfaces. This makes the shell thinner at the top and thicker at the bottom, 
leading to wall non-uniformity. 

♦Appendix IV was authored by R. B. Jacobs of R. B. Jacobs Associates, 
consultant to current contract. 
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We assume that the shell does not rotate about any axis, and that. If It 
has any motion, the motion will be parallel to the uniform acceleration field. 
These assumptions are credible because tumbling of the shell tends to be 
damped. T‘'»n, because of the initial symmetry of the shell, the shell will 
always be axially symmetrical with respect to its diameter that is parallel 
to the acceleration field. It follows that the thickness ( 6 ) of the shell 
is a function of only time (t) and one spatial variable, the angle (J) shown 
in Figure 1 . 

A fundamental assumption in the following development is that the shell 
be thin. That is, if the radius of the exterior spherical surface be R 
and the thickness of the shell be 6 , we assume that 

f <<i 0) 

Note that the shell thickness, 6 (t, 4>), is defined as a length, measured along 
a radius of the exterior spherical surface, between the interior and exterior 
spherical surfaces. It is apparent from the foregoing discussion that the 
minimum shell thickness will be at the top (at ()> = 0) and the maximum shell 
thickness will be at the bottom (at ({> = it). At 41 = 0, 6(t,<j>) 5 6Q(t); at 
4) - ir, 6(t,4>) = non-uniformity (WNU) is defined as: 



In order to proceed with our development we must make decisions regarding 
the conditions at the boundaries of our shell. Two most tenable conditions 
are that the shearing stresses are the same on each side of the interior 
surface, and are the same on each side of the exterior surface. Because the 
liquid shell is bounded by gas, both on the inside and on the outside, it 
would probably be satisfactory to assume that these stresses are zero at 
both surfaces. Even these conditions lead to the unacceptable complexities 
mentioned above. We therefore make assumptions that will yield the simplest 
model, and will also permit the development of wall non- uniformity controlled 
by viscous and acceleration effects: 

1 . There are only tangential velocity components. 

2 . The shearing stress at the interior surface of the shell are zero. 

3 . The velocity of the liquid at the exterior surface, relative to that 





Figure 3 FLUID FLOWS TO CONTROL VOLUME 



surface, is zero (no slip). Conceptually, we imagine that the 
shearing stresses, resulting from the fall of the shell through 
its surrounding medium, are sufficient to accomplish this. 

4. At every instant, a steady-state velocity distribution exists 

across the shell wall; however, this distribution varies with time. 
This assumption is conf.ervative relative to the purposes of this 
development. It leads to predicted values of WNU that are excessive; 
actual shells showed smaller non-uniformities than predicted. 

Applying Newton's Second Law, in the tangential direction, to the fluid 
element in Figure 2 gives 


^ = - p a sin (J) (3) 

where 

T = shearing stress in the tangential direction, 

Z = a radial coordinate, measured inward from the exterior spherical surface, 
p = density of shell material, and 
a = acceleration. 

Assuming that the shell material is Newtonian, 


T = 


dv 
^ dZ’ 


(4) 


where 

u = absolute viscosity, and 
V = velocity in tangential direction. 

Substituting Equation 4 into Equation 3, and utilizing assumptions 2 and 3 in 
the preceding paragraph, we have 

V =^sin (D Z [6 -|]. (5) 

Applying continuity considerations to the control volume in Figure 3 
and noting that 


6 

W((j)) p(2 TT R sin (j> dZ)v, (6) 

Z=0 


we obtain 
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^ fsIft h 


(7) 


We define two dimensionless parameters: 

1. Dimensionless wall thickness, n 

n = (f) (8) 

2. Dimensionless time, x 

t t. (9) 

Introducing Equations 8 and 9 into Equation 7 yields 


-iQL_= . 
3t 


1 

sin (|) 




sin^ (})) 


(10) 


Equation 10 is our final non-linear, partial differential equation which, 
when used with specific initial and boundary conditions, will permit us to 
obtain the shell-wall thickness as a function of time and position: 
i.e., n(x,<f)). 

By utilizing the assumptions made above, most particularly Equation 1 
and the sphericity of the bounding surfaces, Equation 10 can be reduced to an 
ordinary differential equation. Let S(x) be the distance between the centers 
of the interior and exterior spherical surfaces. Then, because n « 1 and 

C 

(•^) f V « 1, we may derive that 


n(x,(|)) 



- v(x) cos (j) L 


where 

= radius of the interior spherical surface. 
Substituting Equation 11 into Equation 10, we obtain 


( 11 ) 


cos <p 


Ri 2 

3(1 - ^ — V cos <p) V sin^ (Ji 

1 

+ 2(1 - V cos (})) * cos (p . 


dv 

dx 


(12) 


Equation 12 is an ordinary differential which may be integrated at 4i » 0 
and (j) » TT. Performing the integrations, and substituting into Equation 11, 
we obtain 
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( 13 ) 


and 



where 

riQ * value of n at (|) » 0, 

n = value of n at (b = IT 

- = value of n when shell wall is uniform, 

n 

Note that 


n = 



(14) 


(15) 


where 

(AR) = aspect ratio of the shell. 

(AR) = (16) 

6 


where 

0 - outside diameter of shell and 


5 * shell wall thickness when shell wall is uniform. 

Substitute Equations 15, 16, and 9 into Equations 13 and 14. Combine 
Equations 8 and 2 to express WNU in terms of tiq and Combining the 
the resulting three equations and rearranging gives: 



(17) 


where 


T * 


8oaD 

MW 


t. 


(18) 


Equations 17 and 18 express the wall non-uniformity (WNU), as a function 
of time (t), of all thin-walled shells which conform with the assumptions 
and limitations delineated in the development. We see that the wall non- 
uniformity is completely specified by the single, dimensionless parameters 
(T) defined by Equation 18. 
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Equation 17 is plotted in Figure 4, where T is considered as a 
dimensionless time. We see, for example, that WNU = 10% when T = .095. 

We note that a characteristic time, t*, can be defined for the system. 


t* 



(19) 


t* is, of course, determined by the shell geometry, shell material properties, 
and the acceleration. Examination of Equations 17, 18, and i9 shows that t* 
is the real time required for the wall non-uniformity to become infinite. 
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